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We review the results of the spin-fermion model for correlated electron materials that are suffi- 
^\ • ciently close to an antiferromagnatic instability that their staggered static magnetic susceptibility 

in the normal state is large compared to that found in a conventional Fermi liquid. We demonstrate 
that for such materials magnetically-mediated superconductivity, brought about by the exchange 
D ' of spin fluctuations, is a viable alternative to conventional phonon-mediated pairing, and leads to 

^ I pairing in the d^2_y2 channel. If the dominant interaction between quasiparticles is of electronic 

origin and, at energies much smaller than the fermionic bandwidth, can be viewed as being due to 
the emission and absorption of a collective, soft spin degree of freedom, the low-energy physics of 
\ these materials is accurately described by the spin-fermion model. The derived dynamic magnetic 

susceptibility and quasiparticle interaction coincide with the the phenomenonological forms used to 
fit NMR experiments and in earlier Eliashberg calculations. In discussing normal state properties, 
, the pairing instability and superconducting properties, we focus our attention on those materials 

^ ' that, like the cuprate, organic, and some heavy electron superconductors, display quasi-two dimen- 

O ' sional behavior. In the absence of superconductivity, at sufficiently low temperatures and energies, 

a nearly antiferromagnetic Fermi liquid is unconventional, in that the characteristic energy above 
which a Landau Fermi liquid description is no longer valid is not the Fermi energy, but is the much 
I smaller spin-fluctuation energy,aJsf. For energies (or temperatures) between uj^f and the Fermi en- 

^ . ergy, the system behavior is quite different from that in a conventional Fermi liquid. Importantly, it 

' is universal in that it is governed by just two input parameters -an effective spin-fermion interaction 

, energy that sets the overall energy scale, and a dimensionless spin-fermion coupling constant that 

diverges at the antiferromagnetic quantum critical point. We discuss the pairing instability cased 
by the spin-fluctuation exchange, and "fingerprints" of a spin mediated pairing that are chiefiy as- 
sociated with the emergence of the resonance peak in the spin response of a d-wave superconductor. 
We identify these fingerprints in spectroscopic experiments on cuprateb superconconductors. We 
conclude with a discussion of open questions associated primarily with the nature of the pseudogap 
. state found in underdoped cuprates. 
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The identification of the microscopic mechanisms re- 
sponsible for superconductivity and the nature of the 
superconducting pairing state continues to represent 
one of the moat exciting theoretical challenges in the- 
oretical physical. In the so-called conventional super- 
conductors, at frequencies less than or comparable to 
the Debye frequency, the attractive phonon-induced in- 
teraction between electrons wins aut over the repul- 
sive screened Coulomb interactioncl and brings about 
superconductivityu. The pairing of electrons in the su- 
perconducting state is in an s— wave channel. The pri- 
macy of phonon-induced interaction in conventional su- 
perconductors has been demonstrated with great clar- 
ity. The phonon density of states, obtained by inelastic 
neutron scattering experiments, and the spectrum of the 
bosons which mediate pairing, as deduced from tunnel- 
ing experiments, agree very "Wj^lt in systems like PbQ'Q. In 
addition to the isotope effectO'Q, this comparison of two 
independent experiments is generally considered to be a 
very reliable proof of a phonon-mediated pairing state. 

The analysis of the tunneling data relies heavily on 
the existence of strong coupling effects in the quasipar- 
ticle density of states and assumes the -ualidity of the 
Eliashberg approach to superconductivityu. Eliashberg 
theory for conventional superconductors is extremely ro- 
bust due to the decoupling of typical electron and phonon 
time scales caused by the small ratio of the velocity of 
sound and the Fermi velocity. This smallness also implies 
that the Debye frequency is much smaller than the Fermi 
energy, and hence the quasiparticles that participate in 
the pairing are low-energy quasiparticles, located in the 
near vicinity of the Fermi surface. Landau Q showed 
that the low energy properties of a normal Fermi liq- 
uid are characterized by a small number of parameters 
and are independent of the details of the underlying lat- 
tice Hamiltonian. Both, the normal and superconducting 
states of a conventional superconductor may be viewed as 
protected states of mattcillj, states whose generic low en- 
ergy properties, insensitive to microscopic details at large 
energy, are determined by a higher organizing principle. 
In this view, in conventional superconductors, the super- 
conducting transition marks a transition from a Landau 
Fermi liquid quantum protectorate to a BCS quantum 
protectorate. The success of the BCS-Eliashberg theory 
for conventional phonon mediated superconductors is al- 
most unique for interacting many body systems. 

It is well known however that the pairing state in the 
Bardeen-Cooper-Schrieffer theory does not need to be 
caused by the interaction between electrons and lattice 
vibrations. Generally there are two distinct classes of the- 
ories of unconventional pairing. The first, and more con- 
servative approach is to replace phonons by another col- 
lective bosonic excitation of the solid. This appmach suc- 
cessfully describes the physics of superfluid "^Hellj, where 
the intermediate bosons are failed ferromagnetic spin 
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fluctuations (ferromagnetic paramagnons) ; a large value 
of the ferromagnetic correlation IsMth is not required in 
view of the Kohn-Luttinger effectiiil. This magnetically- 
mediated interaction causes pairing in a state with an- 
gular momentum I = 1 (p-wave pairing) and leads to a 
rich phaise diagram and a large class of new physical phe- 
nomenaE3~o. Magnetically-mediated superconductivity 
has been proposed for various organic and_heavy fermion 
superconductors by a number of authorstj'Eil. 

A second approach to unconventional pairing is more 
phenomenological, and is based on the assumption that 
the superconducting condensation energy is not deter- 
mined by the attractive interaction mediated by some 
boson but rather by the energy gain due to feedback ef- 
fects associated with pairing. The latter may, in princi- 
ple, occur even for a purely repulsive pairing interaction. 
In general, this approach assumes a non-Fermi-liquid be- 
havior in the normal state. Tm)_examples are the inter- 
layer tunneling model of Refs.t^fO and the mid-infrared 
model of Ref. U4. In both models, the expectation value 
of the Hamiltonian is drastically different in the super- 
conducting and the normal states, and the energy gain 
due to pairing apparently cannot be traced back to some 
boson-mediated attraction. 

The two approaches to unconventional pairing are not 
necessarily in contradiction with each other. First, at 
strong coupling, the pairing interaction mediated by a 
low-energy bosonic mode is highly retarded and is a com- 
plex function of frequency, so it is not straightforward to 
determine whether it is repulsive or attractive. Second 
the pairing obviously changes the form of the fermionic 
self-energy and hence affects the kinetic energy. And 
third, if the bosonic mode is itself made out of fermions, 
then the propagator of this mode also changes when the 
system becomes superconducting. This change affects 
the potential energy of the system. Which of the two 
effects is larger depends on the details of the system be- 
havior, but in any case, there are clear similarities be- 
tween a strong-coupling theory which involves a pairing 
boson, and a scenario based entirely on energy gain due 
to feedback from pairing. 

In what follows we adopt the first approach and investi- 
gate the role of failed antiferromagnetic spin fluctuations 
(antiferromagnetic paramagnons) as a possible cause for 
both anomalous normal state behavior and unconven- 
tional i.e. non s-wave superconductivity. This approach 
is chiefly motivated by the physics of the high temper- 
ature cuprate superconductors which have been shown 
to exhibit both highly anomalous normal state behavior 
and an unconventional pairing ptfLte, with angular mo- 
mentum I = 2 (d-wave pairing)E3~EZI. As may be seen 
in Fig.|l], the materials with the highest Tc are located 
reasonably close to an antiferromagnetic state and have 
been shown in nuclear magnetic resonance and inelastic 
neutron scattering experiments to exhibit significant afc 
tiferromagnetic correlations in the paramagnetic state.L^ 
We will shopiria agreement with the calculations of Mon- 
thoux et a/.EiH23, that in a quasi two-dimensional material 



where those correlations are significant (e.g. a spin cor- 
relation length larger than a lattice constant) the normal 
state behavior is anomalous while for Fermi surface pa- 
rameters appropriate for the cuprates, one always gets a 
dx^_y2 superconducting pairing state. We discuss other 
materials below, following a brief historical overview of 
the developments in the spin-fluctuation approach over 
the last decade. References to earlier works can be found 
in the papers cited below. 




doping x 

FIG. 1. Generic phase diagram of high temperature super- 
conducting cuprates. The thermodynamic phases (antiferro- 
magnetic at low doping and superconducting at higher dop- 
ing) are depicted by the shaded regions. The remaining lines 
are either phase transitions or crossovers, visible in a variety 
of experiments. 

A dx2-y2 pairing state in two dimensions due to the 
exchange of near-antiferromagnetic spin fluctuations was 
found in the|-detailed Hubbard model calculations of 
Bickers et alxB. For parameters believed to be rele- 
vant for cuprates in 1987, the superconducting transi- 
tion temperature was comparatively low (< 40 K) under 
what seemed to be optimal conditions. Furthermore, Tc 
decreased as one increased the planar hole concentration 
from a low level, in contrast to experiment. These results, 
when taken together with the early penetration-depth ex- 
periments that supported an s-wave pairing state, were 
responsible for the fact that the magnetic mechanism and 
dx2_y2 pairing had been abandoned by most of the high 
temperature supeijconductivity community by the end of 
1989 (Bedeh et aM). 

At aboutpLhis time, theoretical groups in TokycO'EHl 
and UrbanatilS independently began developing a semi- 
phenomenological, macroscopic theory of spin-mediated 
pairing. Both groups assumed that the magnetic inter- 
action between the planar quasiparticles was responsible 
for the anomalous normal state properties and found a 
superconducting transition to a (1^2 _y2 pairing state at 
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a significantly higher temperature than thpae-achieved 
using the Hubbard model. Moriya et aLt3c3 used a 
self-consistent renormalization group approach to char- 
acterize the dynamic spin susceptibility. The resulting 
effective magnetic interaction between the planar quasi- 
particles was then used to calculaie. Tc and the normal 
state resistivity. Monthoux et a^E3 did not attempt a 
first-principles calculation of the planar quasiparticle in- 
teraction. Rather, they turned to experiment and used 
quasiparticles whose spectra was determined by fits to 
band structure calculations and angular resolved pho- 
toemission spectroscopy (ARPES) experiments. The ef- 
fective magnetic interaction between these quasiparticles 
was assumed to be proportional to a mean field dynanuc 
spin susceptibility of the form developed by Millis et alS3 
that had been shown to provide an excellent descr^j. 
of NMR experiments on the YBa2Cu307_y syste: 

Both groups followed up their initial weak coupling 
calculaticuis with strong coupling (Eliashberg) calcula- 
tionsLJTa that enabled them to take into account lifetime 
effects brought about by the strong magnetic interaction. 
These calculations showed that dx2-y2 superconductiv- 
ity at-iigh Tc is a robust phenomenon. Monthoux and 
PinesEJ also found in a strong coupling calculation that 
they could obtain an approximately correct magnitude 
and temperature dependence of the planar resistivity of 
optimally doped YBa2Cu307_y using the same coupling 
constant (and the same parameters to characterize the 
quasiparticle and spin spectrum) that had yielded a Tc 
of approximately 90 K. They concluded that they had 
established a "proof of concept" for a nearly antiferro- 
magnetic Fermi liquid (NAFL) description of the anoma- 
lous normal state behavior and a spin fluctuation mech- 
anism for high temperature superconductivity. Referring 
back to Fig. 1, these calculations should apply to the 
right of the Tcr-line, where the normal state is an uncon- 
ventional Fermi liquid in which the characteristic energy 
above which quasiparticles loose their Fermi liquid be- 
havior is of order of the spin fluctuation energy and low 
compared to the ferniionic bandwidth. 

Since their calculations unambiguously predicted a 
d^2_y2 pairing state, Monthoux and Pines challenged the 
experimental community to find unambiguously the sym- 
metry of the pairing state. At that time (1991-1992), only 
NMR Knight shift and Cu spiii4a|tt|ice relaxation rate re- 
sults supported (1^2 _y2 pairingp3c3. However within the 
next year or so, the tide turned dramatically awaii-froni 
s-wave pairing, with AREKSa, penetration depthcj and 
new NMR experimentscZI cj on the oxygen spin-lattice 
relaxation time and the anisotropy of the copper spin- 
lattice relaxation time all supporting a (1^2 _y2 state. The 
decisive experiments were the direct phase-sensitive tests 
of pairing symmetry—carried out by Van Harlingen and 
his group in Urbaad^j as well as by Kirtley, Tsuei, and 
their collaboratorsc3. 

In subsequent work on the spin-fiuctuation mechanism 
a microscop|i% Hamiltonian approach to the problem was 
developedE3H. It was shown that the low-energy physics 



of spin-mediated pairing is fully captured by a model that 
describes the interaction of low energy fermionic quasi- 
particles with their psfp, collective spin excitations (the 
spin-fermion model)E2fEa. In particular, it was demon- 
strated that the phenomenological interaction between 
quasiparticles could be derived in a controllable way, even 
at strong coupling, by expanding cither in the inverse 
number of hot spots in the Brillouin zone (— 8 for the 
physical case), or in the inverse number of fermionic fla- 
vors. We discuss this theory in detail in Section 4. As 
will be seen there, the spin-fermion model contains only 
a small number of parameters. These uniquely deter- 
mine system behavior that is fully universal in the sense 
that it does not depend on the behavior of the underlying 
electronic system at energies comparable to the fermionic 
bandwidth. It is therefore possible to verify its applicabil- 
ity by first using a few experimental results to determine 
these parameters, and then comparing the predictions of 
the resulting parameter-free theory with the larger sub- 
set of experimental results obtained at temperatures and 
frequencies which are much smaller than the fermionic 
bandwidth. A major prediction of the spin fermion model 
is that the upper energy scale for the Fermi liquid behav- 
ior progressively shifts down as the system approaches a 
quantum-critical point at T = 0, and there emerges a 
large intermediate range of frequencies where, on the one 
hand, the system behavior is still a low-energy one and 
universal, and on the other hand, it is quantum-critical 
and not a Fermi liquid. 




FIG. 2. The phase diagram of the layered organic super- 
conductor K-(ET)2Cu[N(CN)2]Cl in the units of temperature 
and pressure (from Ref. p5|). PI refers to a paramagnetic 
insulating regime, M to a metallic regime, AF to an antifer- 
romagnetic regime, and SC to a superconducting regime. In 
the region AF-SC, superconductivity and antiferromagnetism 
co-exist. In the region U-SC, the system is an unconventional 
superconductor 
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Cuprate superconductors are not the only candidate 
materials for spin fluctuation mediated pairing and 
non-Fermi liquid quantum protectorates. A number 
of organic superconductors are anisotropic quasi-two- 
dimensional materials that exhibit many of the anoma- 
lies typical of a system with an unconventional pairing 
state. The phase diagram of a quasi two-dimensional or- 
ganic compound k -BEDT-TTF is shown in Fig. |^. One 
can see that, as in the case of the cuprates, the super- 
conducting phase is found in the prmity of an antifer- 
romagnctic phase. Several groupsHtJ have used spin- 
fluctuation theory to predict the position of nodes of the 
superconducting order parameter of these materials. An 
unconventional order-parapjfiter with nodes of the gap-is 
indeed supported by NMEEjLS, thecmal conductivity^, 
millimeter transmissioncil and STME3 experiments. How- 
ever, the last two experiments seem to come to different 
conclusions as. far as the position of the nodes is con- 
cerned. Ref.EHl also finds nodes, but at a position that 
is not consistent with the prediction of a spin fluctua- 
tion induced-pairing state. Finally, penetcation depth 
experimenta23 and recent specific heat datao appear to 
support a conventional s-wave gap. Given those contra- 
dictory experiments, whether quasi-two dimensional or- 
ganic superconductors exhibit an unconventional pairing 
state is, as of this writing, an open question. 




1 2 
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FIG. 3. The phase diagram of Celna in the units of tem- 
perature and pressure (from Ref. p7[ ). 

Cerium-based heavy electron superconductors repre- 
sent another class of strongly correlated electron su- 
perconductors for which a spin-fluctuation induced in- 
teraction between quasiparticles is a strong candidate 
for the superconducting mechanism. Examples include 



CeCusSiaEa, CePdzSia, CelnjB and the newly discos 
ered 1-1-5 materials CeXIns with X = Co, Rh and IrES 
or mixtures thereof. As may be seen in the phase di- 
agrams of Fig.^ and FigJJ, all these materials are close 
to antiferromagnctism, with superconductivity occurring 
close to the critical pressure or alloy concentration at 
which the magnetic ordering disappears. Moreover, ther- 
mal conductivity measurements of CeCoIns, which be- 
comes superconducting at 2.4K at ambient pressures - 
the highest known value of Tc for a heavy fermion based 
system - strongly support a superconducting gap with 
nodes along the (±7r, ±7r) directions, as found in a dx2_y2 
pairing statellj. Another exciting aspect of these systems 
is that by changing the relative compositions of Ir and 
Rh in 1-1-5 materials CeRhi_xIrxIn5, one can move the 
system from an antiferromagnetic to a superconducting 
state at ambient pressure. 




FIG. 4. The phase diagram of Ce2X Ins with X=Co, Rh, 
and ' 
et alh). 



id Ipfn the units of temperature and doping (from Pagliuso 
alM). 



Another widely studied material in which pair- 
ing is ppsfdhly due to spin fluctuation exchange is. 
Sr2Ru04bl"y, where NMR Knight shift experiments 
and spin-polarized neutron scattering measurementsE^ 
reveal that the spin susceptibility is unchanged upon en- 
tering the superconducting state, consistent with spin- 
triplet superconductivity. 

In summary, the cuprates, the 1-1-5 heavy fermion 
materials and the layered organic superconductors are 
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strongly correlated materials that exhibit unconventional 
normal state and superconducting behavior, while the su- 
perconducting phases are located in the vicinity of mag- 
netic instabilities in their corresponding phase diagrams. 
It is then quite natural to assume that in all three cases 
magnetic interactions play a dominant role in the pairing. 

The presence of antiferromagnetic and superconduct- 
ing regions in the phase diagram raises the question 
of whether antiferromagnetism and superconductivity 
should be treated on equal footing in a spin fluctuation 
approach. If they should, the theoretical analysis would 
be complex. Fortunately, this is not the case, at least as 
long as the characteristic energy scales for the magnetic 
interactions are smaller than the fermionic bandwidth. 
The point is that superconductivity is generally a low- 
energy phenomenon associated with fermions in the near 
vicinity of the Fermi surface. On the other hand, antifer- 
romagnetism originates in fermions with energies compa- 
rable to the bandwidth. Perhaps the easiest way to see 
this is to formally compute the static spin susceptibil- 
ity in the random phase approximation (RPA). An RPA 
analysis yields x"^(q) oc 1 - gcff(q)n(q) where 5cff(q) 
is some effective interaction, and n(q) is the static spin 
polarization operator (a particle-hole bubble with Pauli 
matrices in the vertices). For an antiferromagnetic insta- 
bility we need goff(Q)II(Q) = 1. One can easily make 
sure, by evaluating n(Q) for free fermions, that the mo- 
mentum/frequency integration in the particle-hole bub- 
ble is dominated by the upper energy limit that is the 
fermionic bandwidth. This implies that whether or not 
a system orders antiferromagnetically is primarily deter- 
mined by high-energy fermions that are located far away 
from the Fermi surface, and hence the antiferromagnetic 
correlation length, that measures the proximity of a ma- 
terial to a nearby antiferromagnetic region in the phase 
diagram, should not be calculated but rather be taken 
as an input for any low-energy analysis. We discuss the 
practical meaning of this separation of energies in Section 
4. 

A more subtle but important issue is whether the dy- 
namical part of the spin susceptibility should be con- 
sidered simply as an input for a low-energy model (as 
in the case for phonons), or whether the spin dynamics 
is produced by the same electrons that are responsible 
for the superconductivity and hence needs to be deter- 
mined consistently within the low-energy theory. The 
first issue one has to consider here is whether a one- 
band description is valid, i.e., whether localized elec- 
trons remain quenched near the antiferromagnetic in- 
stability and form a single large Fermi surface together 
with the^cmiduction electrons to which they are strongly 
coupledEZl'Cj, or whether the magnetic instability is ac- 
companied by the un-quenching of local moments. In the 
latter case, the volume of the Fermi surface changes dis- 
continuously at the magnetic transition and could e.g., 
cause a jump in the Hall cocfficientEj. The quenching 
versus un-quenching issue is currently sir^jttbject of inten- 
sive debate in heavy fermion materialal^Ej. In cuprates 



the quenching versus un-quenching issue does not seem 
to play a role; it is widely accepted that the formation 
of Zhang-Rice singletfil gives rise to a single electronic 
degree of freedom. Similarly, in organic materials, the 
charge transport in the metallic and superconducting 
parts of the phase diagram is due to the same miss- 
ing electrons in otherwise closed filled molecular orbital 
states. Whether or not the spin dynamics originates in 
low-energy fermions then reduces to the geometry of a 
single, large Fermi surface. For a Fermi surface with hot 
spots, connected by the wave vector at which the spin 
fluctuation spectrum peaks, the low-energy spin dynam- 
ics is dominated by a process in which a collective spin 
excitation decays into a particle-hole pair. By virtue of 
energy conservation, this process involves fermions with 
frequencies comparable to the frequency of a spin exci- 
tation. Consequently, the spin dynamics is not an input. 
If, however, the Fermi surface does not contain hot spots, 
spin damping is forbidden at low-energies and spin fluc- 
tuations are magnon-like propagating excitations. It is 
easy to show that in the latter situation, the full form of 
the spin propagator comes from particle-hole excitations 
at energies comparable to the bandwidth and therefore 
should be considered as an input for the low-energy the- 
ory. 

In this chapter we consider in detail the scenario in 
which the Fermi surface contains hot spots and the spin 
damping by quasiparticles is allowed. Our approach to 
the cuprates is largely justified by the results of extensive 
ARPES and NMR and neutron measurements that indi- 
cate that the Fermi surface possesses hot spots, and that 
spin excitations are overdamped in the normal state. 

Whether or not spin fluctuations are overdamped is 
also of significant conceptual importance for spin medi- 
ated pairing, since this mechanism requires that quasi- 
particles be strongly coupled to the collective spin exci- 
tation mode. At first glance, the undamped (magnon) 
form of the spin propagator appears more favorable for 
spin-mediated pairing than the overdamped form. In- 
deed, if one assumes that the spin-mediated interaction is 
just proportional to the spin susceptibility, the magnon- 
like form is preferable. By the Goldstone theorem, in 
the antiferromagnetically ordered state, the transverse 
spin susceptibility x(q) (that yields an attraction in the 
dx^_y2 channel) even diverges as q approaches the anti- 
ferromagnetic momentum Q, hence the d— wave attrac- 
tion appears to be the strongest. This reasoning, hojW^ 
ever, is in co j rr cct. Schrieffer and his coUaboratorsEBE^ 
and othersfiJH have shown that the Goldstone modes of 
an ordered antiferromagnet cannot give rise to a strong 
d— wave pairing because the full spin mediated interac- 
tion is the product of the spin susceptibility and the 
square of the fully renormalized coupling constant be- 
tween fermions and magnons. The latter vanishes in the 
ordered SDW state at q = Q and this effect exactly com- 
pensates the divergence of the static susceptibility. The 
vanishing of the effective coupling is a consequence of the 
Adler principle which states that true Goldstone modes 
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always decouple from other excitations in a systemEa. 

SchriefFer later arguecS that the near cancellation be- 
tween the enhancement of the spin susceptibility and 
the reduction of the effective magnon-fermion interac- 
tion persists in the paramagnetic state as long as spin 
fluctuations remain propagating excitations.p-This would 
substantially reduce (although not eliminateEll) the spin- 
mediated d-wave attraction. This argument is however 
inapplicable to overdamped spin fluctuations. These are 
not Goldstone modes although they become gapless at 
the magnetic instability. Goldstone modes appea* only 
in the ordered state at the smallest q — Q valuesEJ. For 
near-gapless, but overdamped spin excitations, the Adler 
principle does not work. After all, the damping itself 
is due to the strong coupling of the collective mode to 
fermions. Consequently, the spin-fermion vertex does not 
vanish at the magnetic transition and hence cannot can- 
cel out the enhancement of the d-wave interaction due to 
the increase of the spin susceptibility near Q. Thus, over- 
damped spin fluctuations are better for spin-mediated 
pairing than magnon-like excitations. 

Another aspect of the fact that spin dynamics is made 
out of low-energy fermions is that the retarded interac- 
tion which causes the pairing changes when fermions ac- 
quire a superconducting gap. This feedback from quasi- 
particle pairing on the form of the pairing interaction dis- 
tinguishes pairing mediated by overdamped spin fluctu- 
ations from conventional phonon induced pairing. In the 
latter the bosonic propagator is an input and is only very 
weakly affected by the opening of the gap in the quasi- 
particle spectrum. We will discuss in detail how feedback 
forces one to go beyond an approach in which one solely 
replaces a phonon by a spin fluctuation, and requires that 
one consistently calculates the spin dynamics at low en- 
ergies. While doing this is a theoretical challenge, the 
approach is appealing since it reduces the number of un- 
known parameters in the problem. In particular, we will 
see that in the superconducting state, the propagator of 
spin fluctuations acquires the same form as for optical 
phonons, but the collective mode that is the analog of 
the phonon frequency is fully determined by the super- 
conducting gap and the normal state spin damping. This 
gives rise to new, unique "fingerprints" of spin mediated 
pairing, whose presence can be checked experimentally. 

What is the role of dimensionality? As noted above, 
many of the candidates for spin-mediated pairing are 
strongly anisotropic, quasi-two dimensional systems. 
This not only holds for the cupratcs, but also for a large 
class of organic superconductors. Also, heavy fermion 
superconductors such as CeCoIns display a considerable 
spatial anisotropy. On the other hand, Celna and to a 
lesser extent CeCu2Si2 do not display appreciable quasi 
two-dimensionality in their electronic properties. The di- 
mensionality of the electronic system is important to the 
spin fluctuation model for both normal state and super- 
conducting behavior. We will see that the dynamics of 
the fermions in the normal state is very differently af- 
fected by antiferromagnetic spin fluctuations in two and 



in three dimensional systems. While in the latter case 
only small (logarithmic) corrections to the ideal Fermi 
gas behavior occur in the vicinity of hot spots, we shall 
see that in 2d systems, the strong interaction between 
fermions and spins gives rise to non- Fermi liquid, dif- 
fusive behavior of low energy fermions as the quantum 
critical point is approached. The importance of dimen- 
sionality for superconductivity has been emphasized by 
Monthoux and LonzarichL3 who have shown that it exerts 
a considerable influence on the superconducting transi- 
tion temperature. They pointed out that in three di- 
mensions one cannot avoid repulsive contributions to the 
pairing interaction in choosing a pairing state with nodes, 
so that the same spin-mediated quasiparticle interaction 
is far less effective in bringing about superconductivity 
in three dimensions than in two. 

Since the non- Fermi-liquid behavior of fermionic quasi- 
particles extends down to progressively lower frequencies 
as one approaches the magnetic transition at T = 0, one 
can inquire whether pairing near this quantum-critical 
point is caused by the fermions at the lowest energies that 
are still coherent, or comes from those at higher energies 
(that are still smaller than the bandwidth) that display 
non-Fermi-liquid behavior. If only coherent fermions are 
involved in the pairing, thcii^ according to McMillan's 
extension of the BCS theorjQ, the resulting supercon- 
ducting transition temperature, Tj^, is comparable to 
the upper energy cutoff of the Fermi liquid regime, and 
thus will be of the order of the spin fluctuation energy. 
This energy vanishes at the critical point, and therefore 
magnetic criticality is unaffected by pairing (see the left 
panel in Fig.|^) . If, however, "non- Fermi-liquid" fermions 
can give rise to a pairing instability, then the onset tem- 
perature of this instability in the particle-particle chan- 
nel, that we will identify with Tcr in the phase diagram 
of Fig.^ generally scales with the upper cutoff energy 
for the q uantum-cr itical behavior and remains finite at 
criticalitytH2lllMll23. In this situation, the quantum crit- 
ical point is necessarily surrounded by a dome beneath 
which pairing correlations cannot be neglected as shown 
in the central panel of Fig||. The critical behavior in- 
side and outside the dome is different, and the "primary" 
critical behavior (which gives rise to pairing) can only be 
detected outside the dome. We will demonstrate below 
that the incoherent pairing temperature saturates at a fi- 
nite value when the magnetic correlation length diverges. 
Furthermore, for parameters relevant to cuprates at low 
doping, this temperature is of order of the magnetic ex- 
change interaction J, i.e., it is not small. 

A related issue is whether the pairing instability at 
Tcr implies the onset of true superconductivity (i.e., 
Tcr = Tc), or whether it marks the onset of pseudogap 
behavior. In the latter scenario, for which we will see 
there is considerable experimental support, fluctuations 
prevent a superfluid density ps from developing a nonzero 
value until one reaches a much smaller Tc, and the paired 
incoherent fermions do not participate in superconduc- 
tivity. The phase between Tcr and Tc would then be a 
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neiK-state of matter, the pseudogap state. Abanov et 
aZ.t23 conjectured that a pseudogap regime is a univer- 
sal feature of the spin fluctuation scenario, as below Tcr, 
quasiparticles that are paired into singlets still remain 
incoherent and cannot carry a supercurrent. True super- 
conductivity is reached only at much smaller Tc where the 
systems recovers coherent, Fermi-liquid behavior (see the 
right panel in Fig||). 




q.c. point 



doping 



FIG. 5. The candidate phase diagrams in the units of tem- 
perature and doping of a one-band electronic system near an 
antiferromagnetic quantum critical point. Left panel - the 
phase diagram for the hypothetical situation when only co- 
herent, Fermi liquid quasiparticles contribute to the pairing 
(the McMillan theory applied to spin fluctuations). The anti- 
ferromagnetic and superconducting regions are completely de- 
coupled. Central panel - the solution of the coupled set of the 
Eliashbfpgiequations for the onset of spin-mediated pairing in- 
stabilityt23. The solution shows that at strong coupling, the 
the pairing instability is predominantly produced by incoher- 
ent fermions, and the instability temperature remains finite 
at ^ = 00. Right panel, the proposed phase diagram based 
on the solution of the Eliashberg equations below the pair- 
ing instability-iand general arguments about superconducting 
fiuctuationgi22 . 



The term pseudogap was introduced by FriedelEl^ 



to 

describe the fact that in the underdoped regime of the 
cuprates, the planar quasiparticles begin to develop a 
gap-like structure well above Tc- This behavior was first 
seen in Knigiit-.shift measurements of the uniform spin 
susceptibilitylll3, and later detected in almost all mea- 
sured properties of underdoped cuprates. At present, 
the physics of the pseudogap phase in the underdoped 
cuprate superconductors is not yet fully understood and 
its origin continues to be an open question. We be- 
lieve that the "magnetic scenario" for the pseudogap pro- 
vides a reasonable explanation, but many details still 
need to be worked out. Some researchers on the other 
hand have suggested that the pseudogap phase emerges 



due to strong fluctuatior 
ducting order parameter. 



Qi-tJue phase of the supercon- 
. Others suggest that the 
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pseudogap phase may actually be a new phase of matter 
with. a . hi dden order parameter associated with bond cur- 
rentsEj^Ell (most recently, this, idea has been explored in 
detail by Chakravarty et alSj). A somewhat more gen- 
eral phenomenological possibility discussed by several re- 
searchers is that there exists an additional quantum crit- 



ical point-Qf yet unknown origin slightly above optimal 
doping£3L3 (a number of experiments suggest that this 
point is at doping concentration x ~ 0.19). The pseu- 
dogap and Fermi liquid phases are assumed to be to the 
left and to the right of this new quantum critical point, 
respectively. 

Our main goal is to discuss in detail the "primary" 
quantum-critical behavior within the magnetic scenario 
and how it gives rise to pairing at Tcr- A detailed theory 
of the pseudogap state of high temperature superconduc- 
tors is beyond the scope of this Chapter. However, in the 
interest of providing a base line against which to compare 
both experiment and future theoretical developments, we 
summarize the predictions of the spin-fermion model for 
the pseudogap in Section 6 and discuss other alternatives. 

To spell out the expected regions of applicability to the 
superconducting cuprates ( in doping and temperature) 
of the spin-fluctuation theory without pseudogap physics 
involved, we return to the candidate generic phase dia- 
gram in Fig.|l|. The two lines, Tcr and T* determine dis- 
tinct regimes of physical behavior. Above Tcr, pseudogap 
physics plays no role; the theory of a nearly antiferro- 
magnetic Fermi liquid (NAFL) presented in this Chapter 
should be applicable for both the normal state and the 
superconducting state. Since Tcr crosses Tc near the opti- 
mal doping concentration, the theory with no pseudogap 
involved is roughly applicable at and above optimal dop- 
ing (from an experimental perspective, optimally doped 
materials do show some pseudogap behavior, but only 
over a very limited temperature regime). For the over- 
doped and nearly optimally doped cuprates the transi- 
tion is then from a nearly antiferromagnetic Fermi liquid 
to a BCS superconductor with dj,2_y2 -pairing symmetry. 
We will argue in Section 6 that there is a great deal of 
experimental evidence that at and above optimal doping 
the normal state is indeed a NAFL, and that the pair- 
ing is of magnetic origin. It is also likely that the theory 
can also be extended into g sp-c.alled "weak pseudogap" 



and r*E3ll2j, but we will not discuss 



regime between Tc-, 
this issue here. 

In Section 2 we introduce and motivate the spin 
fermion model that we use to study spin fluctuation in- 
duced pairing. We discuss the weak coupling approach 
to the pairing problem and the symmetry of the mag- 
netically mediated pairing state. In Section 3 we re- 
view the main results and arguments used to justify 
Eliashberg theory for conventional phonon superconduc- 
tors. In particular, we discuss the physical origin of 
the Migdal theorem that allows a controlled approach 
to phonon-induced pairing. In Section 4 we then ana- 
lyze in detail the strong coupling theory for the spin- 
fermion model. We flrst discuss the normal state prop- 
erties of this model and calculate the low frequency dy- 
namics of quasiparticles and spin fluctuations. We next 
consider spin-fluctuation induced superconductivity. We 
show that for magnetically-mediated superconductivity 
one can again analyze the pairing problem in controlled 
calculations that on the level of the equations involved 
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resemble the Eliash.berg equations for electron-phonon 
superconductivityOEil. We demonstrate that the actual 
physical origin of the applicability of a generalized Eliash- 
berg approach for spin mediated pairing is qualitatively 
different from the phonon case, and is associated with 
the overdamped nature of the spin excitations. We solve 
the resulting equations in certain limits and investigate 
the role of quantum critical pairing. In Section 5 we 
present a general discussion of some of the observable 
fingerprints of spin fluctuation induced superconductiv- 
ity, and in Section 6 wc compare our results with ex- 
periments and discuss to what extent the fingerprints of 
spin mediated pairing have already been seen in opti- 
mally doped cuprate superconductors. In our concluding 
Section 7 we summarize our results and comment on sev- 
eral topics that are of interest for a further understanding 
of spin mediated pairing, including the extent to which 
our theory can be extended to address the physics of the 
pseudogap state in underdoped cuprates. 

II. SPIN-FERMION MODEL 

A. Physical motivation of the spin fermion model 

We first discuss the formal strategy one has to follow to 
derive an effective low-energy model from a microscopic 
Hubbard-type Hamiltonian with a four fermion interac- 
tion: 

^ = E ^k^k,a^k,a (1) 

+ E ^klMMM V'ki,aiV'k2,Q2'^k3,Q3V'k4,Q4 

Here U^^'^^'^^'^^ is the four- fermion interaction, ipj^ ^ is 
the creation operator for fermions with spin a and mo- 
mentum k, and £k is the band-structure dispersion. For 
a one band Hubbard model with local Coulomb interac- 
tion, 

^ki,k2,k3,k4 - f^'5ki+k2-k3-k4 (^) 

^ (f^ai Q4 f^a2a3 ^aiQ3^a2Ct4) • 

In a perturbation theory for Eq.|l] involving U and the 
fermion band width, the contributions from large and 
small fermionic energies are mixed. However, near a 
magnetic instability much of the non-trivial physics is 
associated (at any U) with the system behavior at low 
energies. To single out this low-energy sector, one can 
borrow a strategy from field theory: introduce a charac- 
teristic energy cut off, A, and generate an effective low 
energy model by eliminating all degrees of freedom above 
A in the hope that some of the system properties will be 
universally determined by the low-energy sector and as 
such will not depend sensitively on the actual choice of A. 
Eliminating these high energy degrees of freedom is the 



central theoretical difficulty in the field of strongly cor- 
related electron systems. In our case, except for a poorly 
controlled RPA analysis, there are no known ways to per- 
form such a renormalization procedure. Furthermore, the 
separation between high energies and low-energies can 
be rigorously justified only if the interaction is smaller 
than the fermionic bandwidth. For larger interactions, 
the critical exponents very likely will remain the same, 
but the pre-factors will depend on system properties far 
from the Fermi surface. 

Still, it is possible to assume that the separation of 
scales is possible, i.e., that the interaction is smaller than 
the bandwidth and to pursue the consequences of that as- 
sumption. This is the strategy we adopt. By itself, this 
does not guarantee that there exists a universal physics 
confined to low energies. This we will have to prove. 
This also does not mean that the system is in the weak 
coupling regime, as near the antiferromagnetic transition 
we will find a strong, near-divergent contribution to the 
fermionic self-energy that comes from low frequencies. 
What the separation of scales actually implies (to the 
extent that we find universal, low-energy physics) is that 
Mott physics does not play a major role. In particular, 
in our analysis the Fermi surface in the normal state re- 
mains large, and its volume satisfies Luttinger theorem. 
How well this approximation is satisfied depends on dop- 
ing for a given material and also varies from one mate- 
rial to another. Most of our experimental comparisons 
will be made with the cuprates. In cuprates, the Hub- 
bard U in the effective one-band model for Cu02 unit (a 
charge transfer gap) is estimated to be between 1 and 2 
eV. The bandwidth, measured by ARPES and resonant 
Raman experiments, roughly has the same value. This 
suggests that lattice effects, beyond a universal low en- 
ergy theory, do, indeed, play some role. At half-filling, 
lattice effects are crucial as evidenced by the fact that 
half-filled materials are both Mott insulators and anti- 
ferromagnets with local (nearest-neighbor) spin correla- 
tions. Doping a Mott insulator almost certainly initially 
produces a small Fermi surface (hole or electron pockets) . 
This small Fermi surface evolves as doping increases and 
eventually transforms into a large, "Luttinger" Fermi sur- 
face. How this evolution actually occurs is still a subject 
of debate. From our perspective, it is essential that at 
and above optimal doping, all ARPES data indicate that 
the Fermi surface is large. Correspondingly, magneto- 
oscillation experiments in BEDT-TTF based organic su- 
perconductors also show that the Fermi surface of these 
materials is large. We believe that in this situation, lat- 
tice effects change the system behavior quantitatively but 
not qualitatively, and the neglect of lattice effects is jus- 
tified. We emphasize however that our analysis certainly 
needs to be modified to incorporate Mott physics close 
to half-filling. 

Several aspects of our approach have a close similar- 
ity to the fluctuation exchange approximation (FLEX), 
which in case of a single band Hubbard model corre- 
sponds to a self consistent summation of bubble and 
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ladder diagranisll23ll23. Specifically, the emergence of 
a sharp gapless resonance mode in the spin excitation 
spectrum of a d-wave superconductor, the feedback of 
this mode on the fermions and the anomalous nor- 
mal state behavior of low energy fermions close to an 
instability are very similar in both 



qtie, 



antifcrroms 
approachcal23T£2 



On the other hand, the FLEX ap- 
proach attempts to determine the static spin response 
and thus the actual position of the quantum critical 
point in terms of the bare parameters of the model such 
as the local or additional nonlocal Coulomb repulsions 
as well as the band structure Ck- As discussed above, 
the static spin response, characterized by the correlation 
length ^, strongly depends on the behavior of fermions 
with large energy. Details of the underlying microscopic 
model which are hard to specify uniquely as well as un- 
controlled approximations in the treatment of the high 
energy behavior strongly affect the static spin response 
within the FLEX approach, making it hard to discrim- 
inate model dependent aspects from universal behavior. 
It is this latter aspect which is resolved in our approach 
which concentrates exclusively on the universal low en- 
ergy physics for a given ^. 

What should be the form of the low-energy action? 
Clearly, it should involve fermions which live near the 
Fermi surface. It also should involve collective spin 
bosonic degrees of freedom with momenta near Q, as 
these excitations become gapless at the magnetic transi- 
tion. The most straightforward way to obtain this action 
is to introduce a spin-1 bose field S and decouple the 
four-fermi| On | inte raction using the Hubbard-Stratonovich 
procedurelli^JllZl. This yields 

^ = E ^kV-i.^k.a + E C/ (q) Sq • S_q 

k,Q! q 

+ ^ V'k+q,a'^"/3V'k,/3 • S_q (3) 

k,q,a,/3 

where the aaf3 are Pauli matrices and we assumed that 
the four fermion interaction only makes a contribution in 
the spin channel with momentum transfer q. Integrat- 
ing formally over energies larger than A we obtain the 
effective action in the form (see e.g., 118) 



.A 

+5 / V'I+g,aCra/3V'fc,/3 ' + 0(5"*). 

J k,q 



(4) 



The last term is a symbolic notation for all terms with 
higher powers of S. Fortunately, in dimensions d > 2 
these higher order terms are irrele vant , (marginal for 
d = 2) and can therefore be neglect edIIZllIS. 

The integration over k and q in ^ is over 2 + 1 dimen- 
sional vectors q = (q, iuJm) with Matsubara frequency 
ujm- In explicit form, the integrals read 



q-Q|<A (27r) 



in the boson case, and 

Jk J|k-kp|<A {2tt) 



(5) 



(6) 



in the fermion case. Further, g is the effective coupling 
constant, Gq (fc) is the bare low-energy fermion propa- 
gator, and xo (?) is the bare low-energy collective spin 
boson propagator. As we have emphasized, a controlled 
derivation of g, Go(fc), and Xo('?) is impossible. We there- 
fore will not try to calculate g, etc. Rather, we use the 
fact that antiferromagnetism predominantly comes from 
high-energy fermions and further assume that the inte- 
gration over high energies does not produce singularities 
in both bosonic and fermionic propagators. Then, quite 
generally, Go(fc), and Xo(9) should have Fermi- liquid and 
Ornstein-Zernike forms, respectively 



Go{k) 



zo 



lUJri 



Ek 



Xo (q) 



1/ ^ 



(7) 



(8) 



Here, is the bare value of the spin correlation length, 
and the other notations are self-explanatory. The ac- 
tual ^ generally differs from the bare one because the low 
energy fermions that damp the spin fluctuation modes 
might change as well their static properties. ^ acquires 
an additional temperature dependence due to spin-spin 
interactions. We return to this question later. It is es- 
sential that the bare spin propagator does not contain a 
term linear in iv. The latter will only appear when we 
consider the interaction within the low-energy model. 

We also assume that (i) the momentum dependence 
of the effective coupling g is non-singular and can be 
neglected (recall that we are only interested in a nar- 
row range of bosonic momenta near Q), (ii) the low- 
energy fermionic dispersion can be linearized in k — kj?: 
Ek = Vi? • (k — ki?), and (iii) that the Fermi velocity is 
non-singular near hot spots and to first approximation 
its magnitude can be approximated by a constant. In 
doing this we neglect effects due to a van Hove singular- 
ity in the density of states. Finally, the exact values of 
zq and a are not relevant, as both can be absorbed into 
an effective coupling with dimension of an energy 



g = g^z^a 



(9) 



while vf and ^ will always appear only in the combina- 
tion vf£^~^. 

We see therefore that the input parameters in Eq. |^ are 
the effective coupling energy g, the typical quasiparticle 
energy vp^~^ , and the upper cutoff A. An additional pa- 
rameter is the angle (j)o between the Fermi velocities at 
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the two hot spots separated by Q, but this angle does not 
enter the theory in any significant manner as long as it 
is different from or tt. When hot spots are located near 
(0, tt) and (tt, 0) points, as in optimally doped cuprates, 
(j)o is close to 7r/2, the value we use in what follows. 

As we have emphasized, we will demonstrate that the 
low-energy properties of the model are universal and do 
not depend on A, which then can be set to infinity. Out 
of the two parameters that are left, one can construct a 
doping dependent dimensionless ratio 



A 



9 



(10) 



which will turn out to be the effective dimensionless cou- 
pling constant of the problem (the factor 3/47r is intro- 
duced for further convenience). The fact that A scales 
with ^ immediately implies that close enough to a mag- 
netic transition A > 1, i.e., the system will necessarily be 
in a strong coupling limit. Besides A the only other free 
parameter of the theory is an overall energy scale, i.e., 
g (or alternatively the quasiparticle energy vf^~^). All 
physical quantities that we discuss will be expressed in 
terms of these two parameters only. 

Eqs.^H determine the structure of the perturbation 
theory of the model. The interaction between fermions 
and collective spin excitations yields self energy correc- 
tions to both bosonic and fermionic propagators. We will 
show below that at strong coupling, the fermionic self- 
energy strongly depends on frequency and also displays 
some dependence on the momentum along the Fermi sur- 
face. However, its dependence on the momentum trans- 
verse to the Fermi surface can be neglected together with 
vertex corrections. The fermionic and bosonic propaga- 
tors are then given by the Gor'kov expressions, which for 
generality we present in the superconducting state. 



Gk(jt^) = 
Fk(iw) = - 



iu} + Sk(«t^) + Ek 

((iw + Ek(«w))'-$^(H- 
$k(»^) 



Xq (i^) = 



ae 



l + e (q-Q)'-nQ [tLuY 



(11) 



(12) 



(13) 



Here Sk(*t^) and TlQ^iuj) are fermionic and bosonic self- 
energies (k stands for the component along the Fermi 
surface), and i^k(*^) and ^^{iLu) are the anomalous 
Green's function and the anomalous self-energy, respec- 
tively. In the next sections we compute I]k(iw), nQ(icL') 
and $k(iw). 

One can also motivate the spin fermion model by fol- 
lowing the apprjoach developed by Landau in his theory 
of Fermi liquidsO, i.e., by assuming that the influence of 
the other fermionic quasiparticles on a given quasiparticle 



can be described in terms of a set of molecular fields00. 
In the present case the dominant molecular field is an ex- 
change field produced by the Coulomb interaction U{q). 
However, one has to consider this field as dynamic, not 
static. The corresponding part of the action contains 



S 



(14) 



with fermionic spin density s., = ipl^g^^aafjipk.is- As 
in the Landau theory of Fermi liquids, we assume that 
the molecular field HJ^"* is given by the linear response 
function. 



.9^ XO (q) S-q 



(15) 



This expression is valid as long as one is not in a region 
so close to a magnetic instability that nonlinear magnetic 
effects play an important role. Formally, this expression 
can be obtained from Eq.^ by integrating out the col- 
lective degrees of freedom S. A relation between this 
purely fermionic approach and the bosonic spin suscep- 
tibility X (q) of Eq. ^ can be established by evaluating 
the reducible four point vertex in the spin channel in the 
lowest order of perturbation theory. We find 

^al3,js{k, k', q) = ~Vcff{q) (Ja/SCTjS (16) 

where the effective quasiparticle interaction is propor- 



tional to the renormalized spin propagator of Eq. 13 

g'ae 



Ves{q) = g'xiq) = 



i + e2(q-Q)2-nQ(c^)- 



(17) 



B. Weak coupling approach to the pairing instabihty 

One of the most appealing aspects of the spin- 
fluctuation theory is that it inevitably yields an attrac- 
tion in the dx^_y2 channel. As with any unconventional 
pairing, this attraction is the result of a specific momen- 
tum dependence of the interaction, not of its overall sign 
which for spin fluctuation exchange is positive, i.e., re- 
pulsive in the s-wave channel. This differentiates spin in- 
duced pairing from the pairing mediated by phonons. In 
the latter case, the effective interaction between fermions 
is negative up to a Debye frequency. In configuration 
space the spin- fluctuation interaction, Eq.l7, can easily 
seen to be repulsive at the origin and alternate between 
attraction and repulsion as one goes to nearest neighbors. 
It is always repulsive along the diagonals, and this is why 
a c?2,2_y2 state, in which the nodes of the gap are along 
the diagonals, is the energeticallv-^tfeferred pairing state 
due to spin-fluctuation exchangdliJ. 

Suppose flrst that the spin-fermion coupling is small 
enough such that conventional perturbation theory is 
valid. To second order in the spin-fermion coupling, the 
spin mediated interaction has the following form, Eq. |l^. 
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The total antisymmetry of the interaction imphes that 
the orbital part is symmetric when the spin part is anti- 
symmetric and vice versa. The orbital part of the in- 
teraction at low frequencies has the same sign as the 
propagator of optical phonons at co < lod- However, the 
spin part involves a convolution of the Pauli matrices, 
and is different from phononic SapS-yS — SasSfj-y. Using 

CTa/30'7<5 = ^SasSpj - SapS^S WC find 



(18) 



where Taj,i3S = {Sa/B^^s + SasSf3j)/2 and Saj,i3S = 
{SaisS-yS — Sa5Si3j)/2 are triplet and singlet spin configu- 
rations, respectively. We see that there is an extra minus 
sign in the singlet channel. This obviously implies that in 
distinction to phonons, the isotropic s— wave component 
of the interaction is repulsive, i.e., isotropic s— wave pair- 
ing due to spin fluctuation exchange is impossible. There 
are two other possibilities: unconventional singlet pairing 
for which we will need a partial component of xil) fo be 
negative, or triplet pairing, for which the corresponding 
partial component of xil) should be positive. 



FIG. 6. A representative Brillouin zone and Fermi sur- 
face for a near-optimally doped cuprate superconductor. The 
hot spots are connected by the antiferromagnetic wave vector 
Q = (tt, tt). For a weakly incommensurate magnetic response, 
the number of hot spots doubles but simultaneously the in- 
tensity of the magnetic response at momenta connecting pairs 
of hot spots drops by a factor of two. The net effect then re- 
mains the same as for the commensurate response. The figure 
is taken from |5^. 

It turns out that for a spin susceptibility peaked at 
or near the antiferromagnetic momentum Q, the triplet 
components are all negative, and triplet pairing is im- 
possible. We therefore focus on singlet pairing. The lin- 
earized equation for the pairing vertex F\^{iujm) in the 
singlet channel is 

Fik) = -3g £ F{k')xik~k')Goik')Goi-k'). (19) 

If x(fc — k') — Xk-k' (i^m — i^m) was independent of mo- 
mentum, a solution of this equation would not exist. In 



our situation, however, the interaction predominantly oc- 
curs between fermions that are separated in momentum 
space by Q = (vr,7r) (see Fig.^. 

Since both k and k + Q should be near the Fermi sur- 
face, the pairing predominantly involves fermions near 
hot spots khs and khs + Q. As the hot spots are well 
separated in momentum space, we can change the sign 
in the r.h.s. of Eqjl^ by requiring that the pairing gap 
changes sign between khs and khs + Q- Indeed, substi- 
tuting Fi^^Q^iujm) = —Fi^{iuJm) into Eq.[l9|, we find that 
both sides of this equation have the same sign. The pair- 
ing problem then becomes almost identical to that for 
phonons, and the gap equation has a solution at some 

Altogether, there are N = 8 hot spots in the Brillouin 
zone; these form 4 pairs separated by Q, between which 
the gap should change sign. This still does not specify 
the relative sign of the gap between adjacent hot spots. 
However, if the hot spots are close to (0, tt) and symme- 
try related points, as in optimally doped cuprates, the 
gap is unlikely to change sign between them. This leaves 
as the only possibility a gap that vanishes along diagonal 
directions kx = ±fcy, i.e., it obeys dx^-y^ symmetry. 

We emphasize that although we found that the pairing 
gap necessarily should have dx2_y2 symmetry, its momen- 
tum dependence is generally more complex than simply 
coskx — cos ky. To understand— this , we consider how 
one can generally analyze Eq.^!£3. A standard strategy 
is to expand both the pairing vertex and the interac- 
tion in the eigenfunctions of the representations of the 
D4fi symmetry group of the square lattice. There are 
four singlet representations of labeled Aig, Big, B2g 
and . The basic eigenfunctions in each representation 
are 1 in Aig, coskx ~ cos ky in Big, sin kx sin ky in i?2g 
and {coskx — cos ky) sin kx sin ky in A2g. Other eigen- 
functions in each representation are obtained by com- 
bining the basic eigenfunction with the full set of eigen- 
functions with full D^h symmetry. Obviously, the eigen- 
functions with c?2,2_y2 symmetry belong to Big channel. 
The orthogonal functions in the set can be chosen as 
d„(k) — cos nkx — cos nky. 

One can easily make sure that the pairing problem 
close to Tc decouples between different representations. 
As ^ — !■ cx), the Big components of the spin susceptibil- 
ity diverge while other components remain finite. Obvi- 
ously, the pairing is in the Big channel. At moderate ^, 
all components are generally of the same order, but nu- 
merical calculations show that Big components continue 
to be the largest. We therefore neglect other channels 
and focus only on Big pairing. Still, there are an infi- 
nite number of Big eigenfunctions, and the pairing prob- 
lem does not decouple between them. Indeed, expand- 
ing i^k(wm) and x(k,cj,„) in d„(k) as F^ = J^p^P^pO^), 
x(k — k', iuj—ioj') — Xp(*^~*'^')'^p(k)^p(^') ^^d sub- 
stituting the result into (|l9h we obtainE 



Fn{iuJn 



~3g^T ^ Fp{iuJm') Xn{i(^n 



')Tn.p (20) 
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where 



V 



dPk' d„(kOdp(kO 

4^ KJ^+KF 



(21) 



We see that ah r„^p are nonzero for arbitrary n and p as 
the products dn{k)dp{k) and Go(fc)Go(— fc) = l/((wm)^ + 
(eic)^) are symmetric under D^h. As a result, Eq. EG cou- 
ples together an infinite number of partial components 
Fn- As long as the partial components of the suscepti- 
bility are dominated by momenta relatively far from Q, 
all Xn and hence Fn are of the same order. In some par- 
ticular lattice models, e.g. in the Hubbard model with 
nearest neighbor interaction, U = At and ^ '--^ 2, an RPA 
analysis shows that xi is numerically significantly larger 
than all other x„lZI. In this situation, Fk ~ Fic?i(k), i.e., 
the momentum dependence of the gap should closely re- 
semble the cos kx — cos ky form. In real space, this implies 
that the pairing (even though non-local in time) is local 
and involves fermions from the nearest neighbor sites on 
the lattice. In models with different high-energy physics, 
other Xn may dominate. Clearly, the momentum depen- 
dence of the gap is not universal as long as the pairing is 
local in real space. 

Universality is recovered when the susceptibility be- 
comes strongly peaked at Q. In this limit (in which 
we can also rigorously justify restricting our attention 
to the Big channeltH3), all partial components Xn scale 
as log^ and differ only in sub- leading non-logarithmic 
terms. A straightforward trigonometric exercise shows 
that for near-equal Xn, the pairing gap is very different 
from cos k^ — cos ky-. it is reduced near the nodes and en- 
hanced near the maxima at hot spots with the ratio of the 
gap maximum and the slope near the nodes increasing as 
In this situation the pairing problem is confined to 
hot spots, and the pairing symmetry is inevitably d^2_y2. 



III. SUMMARY OF STRONG COUPLING 
THEORY FOR ELECTRON PHONON PAIRING 

For phonon-mediated superconductivity, Eliashberg 
theory offers a successful approach to study the system 
behavior at strong couphng. It is justified by Migdal's 
theorem which states that the vertex corrections, 6g / g 
and {l/vp)dY.{'k, iLu) /dk, are small due to the smallness 
of the ratio of sound velocity and Fermi velocity stem- 
ming from the smallness of the ratio of the electronic 
and ionic masses. Eliashberg has demonstrated that for 
E(k, zw) « 'Sliuj) and f/tot = 9 + ^9 — 9: the phonon- 
mediated pairing problem can be solved exactly. A re- 
view of Eliashberg t heor y for the electron-phonon prob- 
lem is given in Ref. 11£. As we already pointed out, for 
spin fluctuation induced pairing there is no Migdal the- 
orem simply because the spin propagator is made out of 
electrons and hence a typical spin velocity is of the same 
order as the Fermi velocity. It is therefore natural to in- 
quire whether this implies that an Eliashberg-type treat- 
ment is inapplicable to the spin problem. To properly 



address this issue, we review the case of electron- phonon 
interactions and examine why the smallness of the mass 
ratio affects dE(k, icij)/dk but not dT,{\s., ius) / dui . 

The electron-phonon interaction is generally rather i^j; 
volved and has been discussed in detail in the literature 
For our purposes, it is sufficient to consider the sim- 
plest Frohlich-type model of the electron-phonon inter- 
action in which low-energy electrons are coupled to opti- 
cal phonons by a momentum, frequency and polarization 
independent interaction gcp. Despite its simplicity, this 
model captures the key physics of phonon-mediated pair- 
ing. 

The propagator of an optical phonon has the form 



^(q, ium) 



Wo 



(22) 



Here ujq is a typical phonon frequency, which is of the 
order of the Debye frequency, and Vg is the sound veloc- 
ity. Both ujQ and Vs are input parameters, unrelated to 
fermions. The ratio Vg/vp scales as (rajMYl"^ where M 
is the ionic mass, and m is the electron mass. In prac- 
tice, VgjvF ~ 10"^ |. This is a real physically motivated 
small parameter for the electron- phonon problem. 

Far from structural instabilities, wq ~ Vsa where a is 
the distance between ions. To make the analogy with 
spin fluctuations more transparent, we will consider be- 
low the case when loq <C VsCl, i.e., the correlation length 
for optical phonons is large. This last assumption allows 
us to simplify the calculations but will not change their 
conclusions. 

In the presence of two different velocities, it is not a 
priori obvious what is the best choice of a dimension- 
less expansion parameter in the theory. Two candidate 
dimensionless parameters are 

Acp = g^iATriv.VF)) (23) 
(the factor An is chosen for further convenience), and 



Acp = XcpVs/vp 



i9cp/vFf 



(24) 



For simplicity, we set the lattice constant a = 1. Obvi- 
ously, Acp 3> Acp. In practice, Acp ^ 1 for all reasonable 
gap while Acp can be either small or large. 

Which of the two dimensionless couplings appears in 
perturbation theory? Consider first the lowest-order di- 
agram for the spin-fermion vertex at zero external fre- 
quency and zero bosonic momentum with all fermionic 
momenta at the Fermi surface. Choosing T = and the 
X axis along the direction of v^^^ of an external fermion, 
we obtain the vertex correction in the form: 



yep 



9ep 



9cp 



dijJmd^q 



87r3 J {iuJrn - VpqxY + + (WsQ) 



(25) 



The evaluation of the momentum and frequency integrals 
is straightforward. On performing the integration we ob- 
tain for loq ^ vf, 
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Sg, 



ep 



\ — \ 
-^cp — ^cp 

ffep VF 



(26) 



We see that the vertex correction scales as Acp and is 
small. For ivo ~ Ws, the expression is more complex but 
still is of order A. One can easily check that higher-order 
diagrams yield higher powers of Aepj-i.^., are also small. 
This result, first obtained by Migda£2a, is Migdal's the- 
orem. 

The physics behind this result is transparent and can 
be directly deduced from an analysis of the integrand in 
Eq.^. We see that the product of the two fermionic 
Green's functions yields a double pole at qx — itUm/vF ~ 
ujq/vf- The integral over is then finite only due to the 
presence of the pole in the bosonic propagator which is 
located at a much larger momentum q^ ~ iloq/vs. This 
implies that the electrons that contribute to the vertex 
correction are oscillating not at their own fast frequen- 
cies but at much slower phonon frequencies. At van- 
ishing Vs/vf, fermions can be considered as static ob- 
jects with a typical qx ^ qy ^ ujo/vg. The two dimen- 
sional momentum integration over the product of two 
fermionic propagators then yields 0{l/vp). Simultane- 
ously, the frequency integral over the bosonic propagator 
yields D{q ^ luq/vsjO) ~ 0(1). Combining the two re- 
sults, we obtain Sgcp/gcp — 0(A) as in Eq. 




FIG. 7. A Feynman diagram for the electronic self energy 
due to electron phonon interaction in the Eliashberg theory. 

We next consider the fermionic self-energy. Evaluating 
the diagram in Fig.^ and subtracting the contribution at 
Lu = 0, k = ki?, we obtain 



where 



/(k, iujrr,) = g% 



dflm(Pq 



87r3 uji + ni + {vsqr 
1 1 



(27) 



(28) 



VFqx i{^^7n + f^m) - Cfc - VFqx 



We again have chosen the qx axis along the Fermi velocity 

Vp- 

Suppose for a moment that /(k, iwm) is non-singular 
in the limit k ^ and ajm 0. Then the wave func- 
tion renormalization and the velocity renormalization are 
both expressed via 7(0,0). This quantity can be evalu- 
ated in exactly the same way as the vertex correction. 
Performing the integration, we find 



/(0,0) = -Aep 



(29) 



Substituting this into Eq.^ we find Vp^dT,/dk = 
—dY./d{iuj) = Acp. As Acp ^ 1, this result, if true, would 
imply that both derivatives of the self-energy are small, 
i.e., the system is fully in the weak coupling regime. The 
physical interpretation of this result would parallel that 
for the vertex correction: fast electrons vibrate at slow 
phonon velocities, and this well-out-of-resonance vibra- 
tion cannot substantially affect electronic properties. 

However, this result is incorrect - dS/dw is of order 
Acp, not Aep. The answer lies in the singular behavior of 
/(k, iujrn) in the limit of vanishing k and w„j. To see this, 
we need to compute /(k, ioj) in Eq|2^ at small but finite 
Lu and Ek- An explicit computation yields 



/(k,iw„) = -Aep -t- A 



ILUri 



(30) 



demonstrating that /(0,0) and the hmit of J(k, iw) at 
cj, Ek do not coincide. The additional contribu- 
tion in Eq. |3^ comes from the regularization of the dou- 
ble pole in the integrand in Eq.^ and is in fact math- 
ematically similaL_to the chiral anomaly in quantum 
chromodynamicsll£^ . 

It is clear from Eq,30 that /(k, zwm) has singular low 
frequency and small momentum limits. Substituting 
Eq.pG into Eq.p^, we find that the actual self-energy is 



E(k, IWm) = zAep iOm - \ep{ii^m ~ Ck). 

We see that dE jdiij scales with Aep whereas Vp 



(31) 



^'^dT./dk is 
> 1, this 



proportional to Acp. At strong coupling, Aep 
self-energy gives rise to a strong renormalization of both 
the quasiparticle mass and the quasiparticle wave func- 
tion. Still, vertex corrections and the renormalization of 
the Fermi velocity scale with Aep and are small 

To understand the physical origin of the distinction be- 
tween the frequency and momentum dependence of the 
fermionic self-energy, it is essential to realize that the 
second, O(Aep), term in Eq. |3^ is not caused by real 
electron-phonon scattering. Rather, a careful examina- 
tion of the structure of the expression for /(k, LOra) shows 
that this term accounts for the pole in the fermionic par- 
ticle hole polarization bubble at small momentum and 
frequency. This pole is known to describe a zero sound 
bosonic collective excitation - a vibration of the Fermi 
surface . in .. wh ich it changes its form but preserves its 
volume .11231123 This implies that the fermionic self-energy 
is actually caused by coupling of fermions to their own 
zero-sound collective modes, while phonons just mediate 
this coupling. 

The neglect of vertex corrections and Vp^dT,/dk leads 
to the well known Eliashberg equations for = E(iWm) 
and (in the superconducting state) the pairing self energy 
= $ {iujjn) that also depends only on frequency. For 
phonon-mediated superconductors, these equations have 
the form 



14 



KpT (W„ + En) D {iuJrn-n) 



$2 + (uJn + ^nf 



(32) 
(33) 



This set of equations is solved for a given {icurn) which 
by itself is only very weakly affected by fermions. In 
the normal state and at w ^ tjg, the self-energy Eq.|3^ 
reduces to the first term in Eq.pl|. 

We have reviewed the derivation of these well known 
equations to underline the physical origin of the appli- 
cability of the Eliashberg approach. We see that the 
electron-phonon interaction gives rise to two physically 
distinct classes of interaction processes that contribute 
very differently to the fermionic self-energy and vertex 
corrections. In the first class, fast electrons are forced to 
vibrate at slow phonon frequencies (i.e., phonons are in 
resonance, but electrons are far from resonance). This 
gives rise to both vertex and self-energy corrections, 
but these are small (cx Vg/vp) and are neglected in the 
Eliashberg theory. In the second class, phonons mediate 
an effective coupling between fermions and their zero- 
sound collective modes. This process involves fermions 
with energies near resonance and phonons far away from 
resonance, and yields a strong renormalization of the 
fermionic propagator. The applicability of the Eliash- 
berg approach is then based on the fact that one has to 
include the influence of phonon-mediated scattering on 
zero-sound vibrations, but can neglect the direct scatter- 
ing of electrons by phonons. 

These insights into the origin of the applicability of the 
Eliashberg theory will now be used to justify a general- 
ized Eliashberg approach for spin mediated pairing. 



IV. STRONG COUPLING APPROACH TO 
SPIN-FERMION INTERACTION 

As shown in the previous section, for electron phonon 
interactions the smallness of vertex and velocity renor- 
malizations is caused by the small ratio of the Bose ve- 
locity and the Fermi velocity. The spin problem is quali- 
tatively different. The bare spin-fluctuation propagator, 
Eq.^, describes propagating magnons whose velocity c 
is expected to be of the same order as the Fermi veloc- 
ity. There is then no a priori reason to neglect vertex 
and velocity renormalizations. Fortunately, this argu- 
ment is not correct for the following reasons: First, as 
we just found in case of electron-phonon interaction, the 
fermionic self-energy in the Eliashberg theory is insen- 
sitive to the ratio of sound and Fermi velocities. The 
small ratio oivs/vp is only necessary to eliminate regular 
terms in the self-energy, which are due to real scattering 
by phonons. For these terms to be small it is sufficient 
that bosons are slow modes compared to fermions. Sec- 
ond, the dynamics of the spin fluctuations is drastically 



modified, compared to the ballistic behavior of Eq. ||. A 
strong spin-fermion interaction at low-energies changes 
the bosonic dynamics from propagating to diffusive. Dif- 
fusive modes have a different relationship between typical 
momenta and energies compared to ballistic ones, making 
them slow modes compared to electrons. Consequently, 
regular terms in the fermionic self-energy again become 
smaller than singular ones. 

We will show that in dimensions between d = 2 and 
d = 3, dE (cj) /did scales as A'^"'^ where, we recall, A esc ^ is 
the dimensionless spin-fermion coupling, while Vp^dl^/dk 
remains non-singular at ^ = cxd in d > 2 and only log- 
arithmically increases for d — 2. This implies that an 
Eliashberg-type approach near a magnetic instability is 
fully justified for d > 2 and is "almost" justified for d = 2. 
In the latter case (which is the most interesting because 
of the cuprates) , we will have to invoke an extra approx- 
imation (an extension to large N) to be able to perform 
calculations in a controllable way. 

Our strategy is the following. We first establish that 
one can indeed develop a controlled approach to the spin 
fiuctuation problem in the normal state, see also Ref. 
p^ . Next we apply this theory to the pairing problem 
and show that there is indeed d— wave superconductivity 
caused by antiferromagnetic spin fluctuations. We dis- 
cuss the value of Tc near optimal doping and the prop- 
erties of the superconducting state, particularly the new 
effects associated with the feedback from superconduc- 
tivity on the bosonic propagator, since these distinguish 
between spin- mediated and phonon-mediated pairings. 



A. Normal state behavior 

For our normal state analysis we follow Ref. |l] and 
perform computations assuming that the Eliashberg the- 
ory is valid, analyze the strong coupling results, and then 
show that vertex corrections and w^^dEk (w = 0) /dk are 
relatively small at strong coupling. 

We begin by obtaining the full form of the dynami- 
cal spin susceptibility as it should undergo qualitative 
changes due to interactions with fermions. The self- 
energy for the spin susceptibility (that is the spin po- 
larization operator) is given by the convolution of the 
two fermionic propagators with the momentum difference 
near Q and the Pauli matrices in the vertices. Collecting 
all combinatorial factors, we obtain: 



Ilq {iujn 



~8ge 

X Gk+q 



d'^kdflr. 



(2^)3 ^•^(*"™) 



Here Gk {ii^m) is a full fermionic propagator 

1 



Gk (iuJm) = 



E {iuJr, 



(34) 



(35) 



and E {iLUm) is the self energy which remains to be deter- 
mined. In principle, even in Eliashberg theory, this self- 
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energy depends on the momentum component along the 
Fermi surface. However, for computations of IIq {icum) 
this dependence can be safely neglected; both fermions 
in the fermion polarization bubble should be close to the 
Fermi surface, and hence the momentum integration is 
necessarily confined to a narrow region around hot spots. 

The 2 + 1 dimensional integral over momentum and fre- 
quency in Eq. ^is ultraviolet divergent; this implies that 
its dominant piece comes from highest energies which we 
have chosen to be 0{A). This contribution, however, 
should not be counted as it is already incorporated into 
the bare susceptibility. In addition, IIq (iujm) contains 
a universal piece which remains finite even if we extend 
the momentum and frequency integration to infinity. The 
most straightforward way to obtain this contribution is to 
first integrate over momentum and then over frequency, 
i.e. perform computations assuming an infinite momen- 
tum cutoff but keeping energies still finite. One can easily 
make sure that for this order of limits, the high-energy 
contribution is absent due to a cancellation of infinities, 
and the full IIq (itOrn) coincides with the universal part. 
The integration is straightforward and on performing it 
we obtain 



and Pine to describe the spin dynamics observed in 
NMR experiments on the cuprates. 

We emphasize that the polarization operator, Eq.^, is 
independent of the actual form of the fcrmionic self en- 
ergy Skhs i^)- The explanation of why a k— independent 
fermionic self-energy does not affect the polarization 
bubble ai — a finite momentum transfer was given by 
Kadanofl££j who analyzed a similar problem for phonons. 



He pointed out that the linearization of the fermionic 
energy around kp is equivalent to imposing an approxi- 
mate Migdal sum rule on the spectral function A{k, uj) — 
(I/tt) ImG(k,w) 



J dekA{k,uj) = 1, 



(40) 



where et, the bare fermionic dispersion relation, = Vk • 
(k — kp) in our case. Expressing G(k, in terms of 
A(k, u) via the Kramers-Kronig relation and making use 
of Eq.i^, one finds 

n(q,c.) = -^rd.'^-f0(.^), (41) 



IIq (iujjn) = 



(36) 



where the characteristic energy scale is the spin fiuctua- 
tion frequency. 



4 5 sin 00 



(37) 



We recall that 0o is the angle between the Fermi velocities 
at the two hot spots separated by q ~ Q. In optimally 
doped cuprates, 00 ~ tt/2 and weakly depends on q as 
long as one is far from momenta connecting Fermi points 
along the zone diagonal. This implies that we can safely 
set 00 — 7''/2 and obtain, on making use of Eq.nG, 



4 9 



_3_vfC_ 
16 A 



_9_ 

64^ 



(38) 



An analytical continuation of Eq.^ to the real fre- 
quency axis yields Hq (lo) = i|ti;|/cL'sf. As we anticipated, 
the coupling to low-energy fermions gives rise to a finite 
decay rate for a spin fluctuation. As typical spin frequen- 
cies are of order w_f^~^ ^ c^~^, at strong coupling, A > 1, 
the induced damping term uj/uJsf '■^ A uj/{vp£^^^) is large 
compared to the frequency dependent ^^^/(c^^^)^ term 
in the bare susceptibility, i.e., the interaction with low 
energy fermions changes the form of the spin dynamics 
from a propagating one to a relaxational one. Neglecting 
the bare frequency term, we obtain 



where f{uj') is the Fermi function. Eq.^ then follows 
from the fact that /{lli') is 1 at tj' = — oo and at 

LU' = +O0. 

We next determine the fermionic self energy Sk (iuj). 
In Eliashberg theory it is given by 



/(Pqde 
^^^^3 Xq (ie) Gk+q {iujjn + ie) 



(42) 



Consider first the self-energy at a hot spot. Introducing 
jSn = iujn + Ski,^ (ii^n) and q — Q = q, we obtain from 
Eq.il 
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duJm'Cpq 



{2Tif J q2+^-2(^i 



(43) 



with q^ = g| 4- . The integration over the component 
g|| , parallel to the velocity Vk,,^ at khs , is elementary and 
yields 



3ff 



(44) 



Xq 



Z0+ = 



1 + C^(q-Q) -ioj/ujsi 



(39) 



The same purely relaxational dynamic spin susceptibil- 
ity was introduced phenomenologically by Millis, Monien 



where go = ^"Ml + and q_E = I]„+™//u_f. 

We next assume, and then verify, that a typical q± <^ 
qa. q\ in the first term can then be ignored, and the 
momentum integration yields 
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dq 



1 1 1 /"^ ,~ 1 

— 7^=^= ~ — / "9^7 7 

V?! + 9o ^-L - 90 J -co i<ll- - ^1e) 

= iTTsignqE 

= i7rsign(wm + uj'^). (45) 



On substituting this result into Eq,44 and splitting the 



frequency integral into several pieces, depending on the 



sign of ijJr, 



we obtain 



duj' 



1 



(46) 



where A is given by Eq, 10, The remaining frequency in- 
Ek (iujm) = iA , , ■ (47) 



tegration is elementary and yields 

2lu, 



1 



1 



We now return to the approximation we made above, 
that a typical q± <^ qo- Since the typical values for q±, 
i.e. those which dominate the integral in Eq.^, are of or- 
der qE, the approximation is valid if qo 3> qs- Further, a 
typical internal frequency is of order ujm, hence qs ^ 
\ujm. + S(a;„j)|/wF and qo {uJsfS.'^y^Qi^ml + t^sf)"^/^ ~ 



UJr, 



,1/2 



/vp- At weak coupling, A ^ 1, the 
condition qo 3> qs implies that lo <^ 5/A, i.e., at en- 
ergies smaller that g, the approximation is well justi- 
fied. At strong coupling and lu ^ cjgf, qp ~ Xuj/vp, 
qa ~ XuJsf/vp, and the criterion go ^ is again satis- 
fied. Finally, at strong coupling and uj > lUs/, qE ^ Qoi^ 



xl/2 



jvp. In this limit, the approximation we made is 
qualitatively but not quantitatively correct. In order to 
develop a well contcol]cd-,theoretical framework for this 
limit, Abanov tt a/.tZ3'E£3 developed a controllable 1/A^ 



approach by extending the model to a large number of 
hot spots, N ^ in the Brillouin zone [N = 8 in the physical 
case). This allows an expansion in terms of 1/A^. Alter- 
natively, one can extend the model to a large number of 
fermion flavors, M, and expand in 1/Af. For large iV, 
<lol<lE ~ N"^ , i-e. the approximation go ^ <1e is justi- 
fied. Another appealing feature of this expansion is 
that within it, vertex corrections and the dependence of 
the self-energy of the momentum component transverse 
to the Fermi surface are also small in and can be com- 
puted systematically together with the corrections to the 
frequency dependent part of the self energyEJ. To keep 
our discussion focused on the key results, we will not dis- 
cuss further the details of the IJN approach. Rather, 
we just emphasize that (i) Eq.^ is quantitatively cor- 
rect even if go ~ IE and (ii) numerically the difference 
between approximate and more involved "exjact" results 
for the fermionic Eki,^ (w) is only few percentEil. 

We next analyze the functional form of Eq.^. We see 
that Ski,, (iw) scales with A and at strong coupling ex- 
ceeds the bare iu term in the inverse fermion propagator 
G^^ {iuj) = iuj— Eki,3(ja;). As was the case with phonons, 
it originates in the scattering on zero-sound vibrations of 



the electronic subsystem, while spin-fluctuations mediate 
the interaction between electrons and their zero-sound 
fluctuations. Further, Y^i^^^^{iuj) evolves at the same typi- 
cal energy iOgf as the bosonic self-energy. This intercon- 
nection between bosonic and fermionic propagators is one 
of the key "fingerprints" of the spin-fcrmion model. 




FIG. 8. The frequency dependence of the normal state elec- 
tron self energy due to spin fluctuation-fermion interaction for 
a quasi particle at a hot spot (from Ref. |5l| ). 

In Figjl we show the behavior of the quasiparticle spec- 
tral function A{ei^, lu) at various Ck- We see that the spec- 
tral weight of the quasiparticle peak rapidly decreases as 
Ek becomes larger than Wgf . 

For small frequencies, lu <C ujsf, the spin susceptibility 
can be approximated by its static form. For the fermionic 
self-energy we find after analytical continuation 



Skh, (lu) = \[uj + 



IUJ \lu\ 
4a;sf 



; [uu < cjsf). 



(48) 



We see that the quasiparticle damping term, although 
quadratic in lu as it should be in a Fermi liquid, scales 
inversely with Wgf, not with the Fermi energy as in con- 
ventional metals. As Wsf vanishes at the critical point, the 
width of the Fermi liquid region, where damping is small 
compared to lu, progressively shrinks as ^ increases. The 
quasiparticle renormalization factor dTj-k^^{uj) / duj\uj=o — 
A (X ^ increases as the system approaches the magnetic 
quantum critical point. The quasiparticle z— factor si- 
multaneously decreases as 



1 



1 



1 



l + A 



(49) 



j=0 



and vanishes at criticality. 

At frequencies above tUsU the imaginary part of the 
fermionic self-energy resembles a linear function of lu 
over a substantial frequency range up to about 8iu;sf , and 
then eventually crosses over to 



[iuu) = isignw {i^\lu\) 



7TTl,.,hl/2 



(50) 
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where 



9 

16^' 



(51) 



Thus remains finite as ^ — *■ c». At the critical point 
^sf 0, the self energy displays the non-Fermi-liquid be- 
havior of Eq.|5^ down to w = 0. A plot of the fermionic 
self-energy is presented in Fig.^. The intermediate quasi- 
linear regime is clearly visible. Note also that the devia- 
tions from Fermi liquid behavior starts already at small 
uj ~ Wsf/2. 
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FIG. 9. The normal state spectral function. Note the ab- 
sence of a quasiparticle peak. This is the consequence of 
the proximity to an antiferromagnetic quantum critical point. 
The figure is taken from Ref Bll 

We now consider how well the Eliashberg approxima- 
tion is satisfied, i.e., whether vertex corrections and the 
momentum dependent piece in the fermionic self-energy 
are relatively small. To do this one needs to evaluate 
I]k(w = 0) at k 7^ khs and the vertex correction 5g/g. 
The details of the derivation can be found in Ref. |5|. We 
have 



and 



Sk {io = 0) 



5l 
9 



log A, 



(52) 



(53) 



where Q{4>o) = 20o/7''. For 00 « 7r/2, Q « 1. We see 
that these two corrections depend only logarithmically 
on the coupling and at large A are parametrically small 
compared to the frequency dependence of the self-energy. 
Furthermore, at large N, both 5g/g and Sk contain an 
extra factor of 1/iV, i.e., scale as (l/iV)logA. We see 
therefore that (l/A^)logA is the analog of the second 
coupling constant Agj, for the phonon case. Just as for 
phonons, the applicability of the Eliashberg theory is re- 
lated to the fact that this second coupling constant is 
much smaller than the primary coupling A. In our case, 



this requires that A ;:!> (log A)/A'^. We however emphasize 
that the smallness of the two "couplings" is not the result 
of the smallness of the velocity ratio but the consequence 
of the proximity to a critical point. Also, in distinction 
to phonons, the second coupling still diverges at the crit- 
ical point and therefore corrections to Eliashberg theory 
cannot be neglected close to the antiferromagnetic tran- 
sition. These corrections, however, have been analyzed 
within a renormalization group approach in Ref. 171 51 
and found to be of minor relevance at intermediate cou- 
pling discussed here. 

It is also instructive to explicitly compute the fermionic 
self-energy in the same way as we did in Section 3 for sys- 
tems with electron-phonon interaction and verify that the 
reason for the dominance of the frequency dependence of 
the self energy is common to that for the phonon case. 
To see this, we introduce, by analogy with Eq. 0^ 



(54) 



E(k, io^m) = {iuJrn - Ck+Q )/(k, IWm ) . 

In the present case, /(k, iojm) is singular 
cj„i 0, and /(k = khs.iw™ 0) 7^ /(k khs,0)El. 
Evaluating /(khs, i^m) in the same way as in the phonon 
case, we find at a; ^ Ugf 



/(khsjiWm) = Ireg + A 



ek-i-Q 



(55) 



where Ireg — O(logA). This form is the same as that 
found for phonons (see Eq. 30). The analogy implies 
that the dominant, 0(A), contribution to the fermionic 
self-energy comes from magnetically mediated interac- 
tions between fermions and their zero sound excitations, 
whereas the actual spin-fermion scattering process in 
which fermions at forced to vibrate at typical spin fre- 
quencies yields a smaller O(logA) contribution. 

Finally, away from hot spots but still at the Fermi sur- 
face, the fermionic self-energy is given by the same ex- 
pression, Eq. as at a hot spot, but with a momentum 
dependent coupling constant Ak and energy scale ix'sf(k) 
which obey 



Ak 

Wsf(k) 



A/(l- 

Wsf (1 



(56) 



where 5k — jkp— khs| is the momentum deviation from a 
hot spot along the Fermi surface. We see that the effec- 
tive coupling decreases upon deviation from a hot spot, 
while the upper energy scale for the Fermi liquid behav- 
ior increases. The increase of Wgf, however, is counter- 
balanced by the fact that ujsf oc sincjio, and 00, whicLt-wje 
had set to be « 7r/2, increases away from a hot spotllHil. 

We see from Eq.^ that the width of the region where 
Sk(ti-') is independent of k (i.e., the "size" of a hot 
spot) depends on frequency. At the lowest frequencies, 
u) < a;sf(k), Sk(w) — AkW, and the hot spot physics is 
confined to a region of width which progressively 
shrinks as ^ increases. However, at frequencies above 
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Wsf(k), Sk(w) Ak(wu;sf(k)-^/^ = A(a;Wsf)"^''^ is indepen- 
dent of k. Accordingly, physical processes that happen 
on these scales are isotropic (apart from the dependence 
on 4)q). In this sense the whole Fermi surface acts as one 
big "hot spot" . 

One can easily perform the above analysis in dimen- 
sions larger than d — 2. One finds that the quasiparticle 
spectral weight behaves for large A and the lowest ener- 
gies to < LL>si as 



(57) 



for 2 < d < 3, and vanishes logarithmically for d ~ 3. 
Correspondingly, in the quantum critical regime we find 
for the self energy along "hot lines" 



(iLu) cx iuj \uj\ 



(58) 



an expression that transforms into Sk^^ (iw) oc 
iwlog(|u;|) for d = 3. Above d = 3 no non- Fermi liquid 
effects result from the proximity to the quantum critical 
point. This demonstrates that many of the effects caused 
by the incoherent nature of the fcrmions are peculiar to 
2d and are considerably less pronounced in three dimen- 
sional systems. 



B. The 



pairing instability temperature 



We now consider the development of the pairin g ins ta- 
bility in the spin-fermion model. We follow Ref. 122. It 



is customary in an analysis of the pairing problem to 
introduce an infinitesimally small particle-particle ver- 



$|^°''(a;) and study its renormaliza- 



tex$(o)_Jc.,-u;) 
tion by the pairing interaction. The corresponding di- 
agrams are presented in Fig. The temperature at 
which the renormalized vertex diverges, i.e., when the 
equation for the full $k(w) has a nontrivial solution at 
vanishing ^^\lu), marks the onset of pairing. 





FIG 
vertex 



JO. Diagrammatic representation for the pairing 
The solid and wavy lines are fermionic and spin 



fluctuation propagators, respectively. 

As noted above, the spin-mediated pairing interaction 
gives rise to dj.2_y2 superconductivity. We argued above 
that near the magnetic instability, the ga p is maximum 
near the hot spots. One can check (see |122| ) that the 
pairing problem is confined to the hot regions in the sense 
that the momentum integration never extends to |k— khs| 
where the momentum dependence of the self-energy or in- 
finitesimal pairing vertex becomes relevant. We can then 
assume that the pairing vertex is flat near the hot spots. 



The underlying d— wave symmetry then implies that the 
gap changes sign between two hot regions separated by 
the antiferromagnetic Q. A separation into hot and cold 
regions is indeed valid only if typical |k — khs| <kp. We 
will see below that the d— wave pairing problem is con- 
fined to frequencies of order g. For these frequencies, the 
width of the hot region is finite and is constrained only 
by the requirement that (k) < g, a, condition that im- 
plies (after a more accurate account of the overall factors) 
that the effective coupling is smaller than the fermionic 
bandwidth. As discussed in the Introduction, the latter 
is a necessary condition for the separation between high 
and low energies, and we assume it to hold. We comment 
below on what happens if typical |k— khs| > kp, i.e., hot 
and cold regions cannot be separated. 

The value of the transition temperature depends sen- 
sitively on the behavior of fermions that are paired by 
the spin-mediated interaction. Our analysis of the nor- 
mal state has shown that the character of the fermionic 
degrees of freedom changes at energies of order Wgf . For 
energies smaller than Wsf, fermions display Fermi liquid 
behavior, while at higher energies they display behavior 
that is different from that in a Fermi liquid. In the BCS 
theory of superconductivity only Fermi liquid degrees of 
freedom contribute to the pairing. Let us suppose that 
this also holds in the present case. Then the pairing prob- 
lem would be qualitatively similar to that of BCS, since 
for frequencies smaller than ajgf, the spin susceptibility 
that mediates pairing can be approximated by its static 
value. The linearized equation for the pairing vertex then 
has the form 



1 + A JTf"L ^ 



(59) 



where the 1 -f A factor in the denominator is the re- 
sult of mass renormalization in the Fermi liquid regime 
(E(a;) ~ Aw). The solution of this equationES yields 



FL 



cjsf exp 



1 + A 
A 



(60) 



At weak coupling, this is just the BCS result. At strong 
coupling, the mass renormalization compensates the cou- 
pling constant, and T^^ saturates at T^^ ~ Wsf . This re- 
sult, if correct, would imply that the pairing fluctuations 
become progressively less relevant as one approaches the 
quantum critical point — > (see the left panel in 
Fig.^). At a first glance, this is what happens, because 
pairing of non Fermi- liquid degrees of freedom seems hard 
to accomplish. Indeed, at frequencies larger than Wgf , the 
pairing interaction decreases as (l-|-|w|/wsf )-i/^ and this 
apparently makes the frequency integral ultraviolet con- 
vergent, i.e., the "logarithmic" pairing problem does not 
appear to extend above Wgf- The flaw in this argument 
is that when the interaction decreases, the mass renor- 
malization produced by the same interaction also de- 
creases, and the large overall A is no longer compensated 
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by 1 + T,{ijj)/uj. Indeed, for u > Wsf, = {ojuj)^/"^, 

where, we recall, u) = (9/167r)g, and the mass renormal- 
ization is 1 + ((I'/tj)^/^ = 1 + 2A(tJsf/w)^/^ < A. Fur- 
thermore, we see that at frequencies between ujsf and w, 
the effective mass and effective interaction both scale as 
The product of the two then scales as I/lj, i.e. 
the "logarithmic" pairing problem extends to frequen- 
cies of the order of Co which, we recall, remains finite at 
^ = oo. 

By itself, this effect does not guarantee that the pairing 
instability temperature is of order uj as the pairing inter- 
action depends on the transferred frequency lo — uj' , and 
the linearized equation for the pairing vertex becomes an 
integral equation in frequency. In particular, for ^ = oo, 
and hence tJsf = 0, we need to solve 



1 



1 



[uj') ■ 

\ ^' V 1^^ ~ ^^'1 



(61) 



Observe that this equation does not have any adjustable 
parameter and is therefore fully universal when T is ex- 
pressed in units of Co. Eq. ^ has been asalyzed in detail 
by Finkel'stein, Abanov and one of ualHS. They found 
that it does have a nontrivial solution at 



to the quasi-two dimensionality. TheJAtter are expected 
to yield Kosterlitz-Thouless physicsE22l. Third, Eq.|6^, is 
only valid when g < W where we recall W ^ vp '^s the 
fermionic bandwidth. In the opposite limit g > W, lat- 
tice effects are important and controlled analytical calcu- 
lations are difficult to perform. One can, however, easily 
estimate that in this limit Tcr oc Vp/g (x lOs/^'^. This es- 
timate coincides with the result of Monthoux and Pines 
who extmcted Tcr oc tOsfS,"^ ^ 'Vp/g from their numerical 
analysiaS. Since is proportional to the hopping ma- 
trix element, and g scales with the Hubbard U at large 
^, it follows that Tcr oc J, where J is the magnetic ex- 
change integral of the corresponding Heisenberg model 
which describes antiferromagnetism at half-filling. 

Eq. |6^ demonstrates that the c?-wave pairing instability 
temperature of a two dimensional system at an antiferro- 
magnetic quantum critical point is finite. It is interesting 
to note that the same holds for p-wave pairing at a fer- 
romagnetic-quantum critical point as shown by Roussev 
and MiUifM. 



C. Superconducting state 



T 



Q.2uj. 



(62) 



They also analyzed the pairing problem at a finite A and 
found that incoherent fermions dominate pairing down 
to a surprisingly small A ~ 0.5. The McMillan like for- 
mula, Eq,60, becomes valid only at smaller A. A plot of 
Tcr versus A is presented in Fig.|ll]. 
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FIG. 11. The results for the instability temperature Tcr ob- 
tained from the solution of the linearized Eliashberg equations 
for different val ues o f the coupling constant A. The figure is 
taken from Ref. 122. 



Several comments are in order here. First, at these val- 
ues of the coupling constant (and doping) Tcr does not 
coincide with the onset temperature for superconductiv- 
ity, Tc, but rather represents the onset temperature for 
pseudogap behavior; the actual Tc is lower, as discussed 
below. Second, we have neglected fluctuation effects due 



We next extend the Eliashberg theory to the spin- fluc- 
tuation induced supercondjicting state. The discussion 
in this section follows Ref.E£3. We derive a generalized 
set of Eliashberg equations for the fermionic self-energy 
and the gap function that include an additional coupled 
equation for the spin polarization operator. The latter, as 
discussed in the Introduction, is produced by low-energy 
fermions and has to be determined self consistently. 

As discussed above, the infinitesimal pairing vertex 
and the fermionic self-energy in the normal state depend 
weakly on momentum in the hot region |k — k^^l < g/W 
where VF oc vj? is the fermionic bandwidth. We will see 
that in the superconducting state, the momentum inte- 
gration is also confined to hot regions. We can then safely 
neglect the weak momentum dependence of both S(ia;) 
and ^{iio), as we did above in calculating Tcr- Subtle 
effects due to this weak momentum dependence will be 
considered in the next section. We will not attempt to 
discuss the behavior of the gap near the nodes. The latter 
is central for the interpretation of the experimental data 
at the lowest temperatures and frequencies, but not at 
energies comparable or larger than the maximum pairing 
gap. From our perspective it is not essential for the spin 
fiuctuation induced pairing state. 



1. Generalized Eliashberg equations 



The derivation of the Eliashberg equations is straight- 
forward. In the superconducting state, the normal 
and anomalous fermionic Green's functions Gk {iujn) and 
Tk {iuJn) and the dynamical spin susceptibility are given 
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by Eqs.|ri| - It is convenient to rewrite Gk and 
Fk (iw„) as 



Gk (iujn) 



Fk (iujn) = i 



(63) 



where il]„ = icij„ + S]kijg {i^n) (in real frequencies, S(a;) = 
w + Ekhs(^))- Without losing generality we can set 
Ek = Vxkx + Vyky and ek+q = -Vxkx + Vyky where 
k = k — khs- The sign change between _Fk and -Fk+q 
is the result of (1^2 -y2 symmetry. The spin susceptibility, 
we recall, is given by 



i + e (q-Q)"-nq {lojm) 



(64) 



We need to obtain the equations for the fermionic self- 
energy Ski,g(«ix'm), the anomalous vertex $kh3(«i^m), and 
the spin polarization operator HQ{iujm)- The spin polar- 
ization operator is obtained in the same way as in the 
normal state, but now there are two particle-hole bub- 
bles: one is the convolution of GkGk+q and the other is 
the convolution of Fk^k+q. We have 

——2 [Gk+q (jw„+™) (65) 
(27r) 

X Gk {iujyn) - Fk+q («w„+,„) Fk {iujni)] 

(the negative sign between the two terms originates in the 
summation over the spin components). The momentum 
integration can be performed explicitly and yields 



Ilq {iujn) 



F ^ [1 - 5 {iUJra) 9 (iuJn+Tn) 



-/ («t^m) / (j^n+rn)] 



where 



and 



/ («w„ 



I]^-t-a>2(zu;„) 
E^ + <i>2(z^„0 



(66) 



(67) 



(68) 



The first term in Eq.|6^ is the result of the regularization 
of the ultraviolet singularity. The additional sign change 
between gg and // terms in Eq|6^ is due to the c?— wave 
form of Fk . 



Eq.^ takes into account the change of the low en- 
ergy spin dynamics in the pairing state. In the case of 
electron-phonon interaction a corresponding change of 
the phonon dynamics exists as well, causing a shift of 
the phonon frequency and line width below Tc. While 
this effect is only a minor corijcetion to the phononic dy- 
namics and is often neglectedllHJ, it leads to a dramatic 
change of the spin dynamics in our case. 

The other two equations are formally the same as for 
phonon-mediated superconductors. The fermionic self- 
energy S(tt'„) is given by 



X (q, ii^-m) Gkh3+q (w„+m) 

(69) 



Performing the momentum integration along the same 
lines as in the normal state calculations we find 

Y^iiUn) ^'^^T'S^ D{iu}ra)g{ii^n+it^7n) (70) 
ZVp ^ — ' 

m 

where D is the effective bosonic propagator that is ob- 
tained by integrating the dynamical spin susceptibility 
over the momentum component along the Fermi surface 
and setting other momentum components to Q (the last 
step is equivalent to the approximation we discussed be- 
low Eq.^). We have 



/dqn 
—X (q, iuJm) 



q± = 



2^1 -Hq (z 



(71) 



An analogous equation is obtained for the anomalous ver- 
tex 

3.9^ 



<I>(iw„) 



2VF 



(72) 



Eqs. 66, [70, and [73 constitute the full set of Ehashberg 



equations for the spin-mediated superconductivity. Al- 
ternatively to $(a;) and T,{uj) we can also introduce 



/ s S(w) , , , $(cj) 



(73) 



The complex function A{uj) reduces to the superconduct- 
ing gap A in BCS theory where we also have Z{oj) = 1. 
In Eliashberg theory, the superconducting gap, defined 
as a frequency where the density of states has a peak, is 
the solution of A{u! ^ A) = uj. 

We again emphasize that the Eliashberg equations are 
valid for fermionic momenta which deviate from hot spots 
by less than g/vp. For these momenta, the pairing ver- 
tex can be approximated by a k-independent function 
which however changes sign between two hot regions sep- 
arated by Q. For larger deviations, the anomalous vertex 
rapidly goes down and eventually vanishes along zone di- 
agonals. 
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2.. Solution of the Eliashherg equations: 

We discuss the general structure of the solutions of 
the set of Eliashberg equations, and then present the 
results of their numerical solution. First, we see from 
Eq. ^ that, as in the normal state, n(a; = 0) = for 
any S(ci;) and This physically implies that the 

development of the gap does not change the magnetic 
correlation length. This result becomes evident if one 
notices that d-wave pairing involves fermions from op- 
posite sub-lattices. Second, the opening of the super- 
conducting gap changes the low frequency spin dynam- 
ics. Now quasiparticles near hot spots are gapped, and a 
spin fluctuation can decay into a particle-hole pair only 
when it can pull two particles out of the condensate of 
Cooper pairs. This implies that the decay into particle- 
hole excitations is only possible if the external frequency 
is larger than 2A. At s niE^ ^ ler , frequencies, we should have 



This result indeed read- 
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n"(w) = at T = 0E£3 
ily follows from Eq. The Kramers-Kronig relation 
n'(w) = (2/7r) J^U"{x)/{x^ - cj2) then implies that 
because of a drop in 11" (w), the spin polarization oper- 
ator in a superconductor acquires a real part, which at 
low bj is quadratic in frequency and is of order oj^/ (Awgf). 
Substituting this result into Eq.|l^, we find that at low 
energies, spin excitations in a d-wave superconductor are 
propagating, gapped magnon-like quasiparticles 



A? 



A2 + c2(q-Q)2-c.2- 



where 



and 



(Acjs 



xl/2 



(74) 



(75) 



' VpA/g. The re-emergence of propagating spin 
dynamics implies that the dynamical spin susceptibility 
acquires a resonance peak which at q = Q is located at 
uj = Ag. 

Eq.^ is indeed meaningful only if A^ < A, i.e., 
Wsf < A. Otherwise the use of the quadratic form for 
n(a;) is not justified. To find out how A depends on the 
coupling constant, one needs to carefully analyze the, .full 
set of Eqs.|66|-|72|. This analysis is rather involvedE£3E"^ 



and is not directly related to the goal of this Chapter. We 
skip the details and quote the result. It turns out that 
at strong couphng, A > 1, i.e. for optimally and under- 
doped cuprates, the condition A > LUgt is satisfied: the 
gap scales with uj and saturates at A « O.SScj = Q.OGg at 
A — > oo, while ujsf oc A^^ 0. In this situation, the spin 
excitations in a superconductor are propagating, particle- 
like modes with a gap A^. However, in distinction to 
phonons, these propagating magnons get their identity 
from a strong coupling feedback effect in the supercon- 
ducting state. 

At weak coupling, the superconducting problem is of 
the BCS type, and A <^ ojgf. This result is intuitively 
obvious as Wsf plays the role of the Debye frequency in 



the sense that the bosonic mode that mediates pairing de- 
creases at frequencies above Wsf. For A <^ cOsf, xiQi^) 
does not have a pole at frequencies where n(a;) oc oj'^. 
Still, .-a-jpoln in . x' \{Q,^) does exist even at weak cou- 
plingtHZlll23^tZ2lll^. To see this, note that at « 2A one 
can simultaneously set both fermionic frequencies in the 
bubble to be close to A, and make both propagators sin- 
gular due to the vanishing of — (f>2 where, we recall, 
t,^uj + S(w). Substituting - ^^(w) cx w - A into 

Eq. ^ and using the spectral representation, we obtain 
for w = 2A -f e 



dx 



{x{e-x)y/^- 



(76) 



Evaluating the integral, we find that H" undergoes a fi- 
nite jump at w = 2A. By the Kramers-Kronig relation, 
this jump gives rise to a logarithmic singularity in TI'^lo) 
at w = 2A: 



, n"(a;) ^, 2A 
ax — cx A log J 



J2A - ^'^ 



2A| 



(77) 



The behavior of n'(w) and n"(a;) is schematically shown 
in Fig.|lJ. The fact that n'(a;) diverges logarithmi- 
cally at 2 A implies that no matter how small A/wgf is, 
x(Q, w) has a pole at Ag < 2A, when n"(a;) is still zero. 
Simple estimates show that for weak coupling, where 
tJsf ^ A, the singularity occurs at As = 2A(1 — Z^) 
where (x e"'^'''/^^^) is also the spectral weight of the 
resonance peak in this limit. 




FIG. 12. Schematic behavior of the real (dashed line) and 
imaginary (solid line) parts of the particle hole bubble in the 
superconducting state. Due to the discontinuous behavior of 
n" (lu) at a; = 2A, the real part is logarithmically di- 

vergent at 2A. For small co, t he real part behaves like uj^/A. 
The figure is taken from 149 



We see therefore that the resonance in the spin suscep- 
tibility exists both at weak and at strong coupling. At 
strong coupling, the resonance frequency is Ag ~ A/A ^ 
A, i.e., the resonance occurs in the frequency range where 
spin excitations behave as propagating magnon-like ex- 
citations. At weak coupling, the resonance occurs very 
near 2 A due to the logarithmic singularity in II'(lu). In 
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practice, however, the resonance at weak coupling can 
hardly be observed because the residue of the peak in 
the spin susceptibility is exponentially small. 

Fig.|l^ shows the results for x(Q: ^) obtained from the 
full solution of the set of three coupled equations at T w 
and three different coupling constantgll£3. For A > 1, the 
spin susceptibility has a sharp peak at w = Ag . The peak 
gets sharper when it moves away from 2A. At the same 
time, for A = 0.5, corresponding to weak coupling, the 
peak is very weak and is washed out by a small thermal 
damping. In this case, x" only displays a discontinuity 
at 2A. 
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FIG. 13. Imaginary part of the dynamical spin susceptibil- 
ity in the superconducting state at T ^ Tc obtained from the 
solution of the set of three Eliashberg equations for coupling 
cons tants A = 0.5, A = 1, and A = 2. The figure is taken from 
149. 



We next show that the resoH|ance peak does not exist 
for s— wave superconductorsEHl. In the latter case, the 
spin polarization operator is given by almost the same 
expression as in Eq. but with a different sign of the 
//-term; recall that the original sign in Eq.(pq) origi- 
nated from the fact that the two fermions in the spin 
polarization bubble differ in momentum by Q, and the 
d-wave gap changes sign under k — ^ k + Q. One can im- 
mediately check that for a different sign of the anomalous 
term, 11" is continuous at 2A. Accordingly, n'(a;) does 
not diverge at 2A, and hence there is no resonance at 
weak coupling. Still, however, one could expect the reso- 
nance at strong coupling as at small frequencies n'(w) is 
quadratic in by virtue of the existence of the threshold 
for n". It turns out, however, that for s-wave pairing the 
resonance is precluded by the fact that n(a-' = 0) changes 
between the normal and the superconducting states. 

One can make sure that in an s— wave superconductor, 
n(u; = 0) < 0, and that this negative term overshadows 
the positive uJ^ term in n(a;) in such a way that for all 
frequencies below 2A, n(a;) < and hence the resonance 
simply does not exist. That n(a; = 0) < in s— wave su- 
perconductors can be easily explained: a negative 11(0) 
implies that the spin correlation length decreases as the 
system becomes superconducting. This is exactly what 



one should expect as s-wave pairing involves fermions 
both from different sub-lattices as well as from the same 
sub-lattice. The pairing of fermions from the same sub- 
lattice into a spin-singlet state obviously reduces the an- 
tiferromagnetic correlation length. 

A.=2, T=0 
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points. This effect accoiuj^ 
of the resonance peak 
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FIG. 14. The real and imaginary parts of the fermionic 
self-energy S(aj) and the pairing vertex <^{uj ) for A = 2 and 
the lowest T. The results are from Ref. [132. 

We also comment on the dispersion of the resonance 
peak. In Eq. ^ we assumed that A^ is a constant. In 
fact, As depends on q since for any given q, A^ cx A(q) 
where A(q) is a d— wave gap at the points at the Fermi 
surface which are connected by q. In particular, A^ 
should vanish at q = Qmin which connects the nodal 

br the "negative" dispersion 
The latter certainly over- 
shadows the positive dispersion due to (q — Q)^ term for 
q close to Qmin and may do so even for q near Q if the 
correlation length is not large enough. This effect is, how- 
ever, not a part of the quantum-critical description (it 
should become progressively less relevant for q ^ Qmiji 
when ^ increases), and we ignore it in the subsequent 
analysis. Note, however, that the negative dispersion of 
the peak implies that the peak exists only for a small 
range of momenta between Q and Qj||i|i„, , In optimally 
doped cuprates, Qmin ~ (O.Stt, 0.87r)PHF^i, and the mo- 
mentum range for the peak does not exceed 6% of the 
Brillouin zone. The actual q region where the peak is ob- 
servable is even smaller as the intensity of the peak also 
decreases when q approaches Qmin- The smallness of the 
q— range for the peak accounts for small overall spectral 
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intensity Iq = J S{c\,uj)(Pqdu> / (27^)^ that turns out to be 
substantially smaller than S{S + l)/?> = 1/4. Still, at Q, 
the intensity of the peak is not small (experirnen-tally, 
] S{(^,Lo)(kj - 1.5 in optimally doped YBCOEBH), 
and wc have verified that for the frequencies that we 
consider below, the typical q — Q that account for the 
feedback on the fermions are well within the q range be- 
tween Q and Qmin- In other words, the small overall 
intensity of the resonance peak does not preclude strong 
feedback effects from the resonance peak on fermionic 
variables. 

X=2., T=0 
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FIG. 15. The real and imaginary parts of the effective gap 
A(tj) and the quasiparticle renormalization factor Z{ijj) for 
A = 2 and the lowest T. The results are from Ref. 
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For completeness, in Figs.^ and ^ we present results 
for the fermionic self-energy and the pairing vertex for 
the smallest T. We see that the real parts of $(w) and 
A(aj) are finite at w = as should be the case in the 
superconducting state. The imaginary parts of $(a;) and 
Yi{lo) (and of A(ijj) and Z{lu)) vanish at small frequen- 
cies and appear only above the threshold frequency that 
is precisely A -|- Ag. Furthermore, all variables have a 
complex internal structure at large frequencies. In the 
next section we discuss the physical origin of the thresh- 
old at A -|- As and also show that one can extract 3A 
from the derivative of I]"(cj). 

Few words about the numbers. For A = 2, A « O.Sw 
and As « 0.2(1;, i.e., A and A^ are comparable to 
each other. For A ^ 1 a numerical solution of the 
Eliashberg equations leads to A '--^ 2Tcr ^ 0.35(D, and 
As ~ 0.25w/A < A. 



V. FINGERPRINTS OF SPIN FERMION 
PAIRING 

In this Section, we discuss the extent to which the 
"fingerprints" of spin-mediated pairing can be extracted 
from experiments on materials that are candidates for a 
magnetically-mediated superconductivity. Due to strong 
spin-fermion coupling, there is unusually strong feedback 
from spin excitations on fermions, specific to d— wave su- 
perconductors with a magnetic pairing interaction. The 
origin of this feedback is the emergence of a propagat- 
ing collective spin bosonic mode below Tcr- This mode 
is present for any coupling strength, and its gap As is 
smaller than the minimum energy ~ 2A that is nec- 
essary to break a Cooper pair. In the vicinity of the 
antiferromagnetic phase. As cx where ^ is the mag- 
netic correlation length. We show that this propagating 
spin mode changes the onset frequency for single par- 
ticle scattering, gives rise to the "peak-dip-hump" fea- 
tures in the quasiparticle spectral function, the "dip- 
peak" features in tunneling SIS and SIN conductances, 
and to singularities and fine structures in the optical 
conductivity. In section 6, we apply these results to 
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fea- 

, . (ii) 
113, tunnelingl^'t^, and conductivity 

are consistent with each other, and (iii) the 
value of As extracted from these various experiments co- 
incides with the resonance frequeiicjf -measured directly 



in neutron scattering experiments 
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A. The physical origin of the effect 

The physical effect that accounts for dips and humps 
in the density of states and spectral function of cuprates 
is not new and is known for convMitioa^l s— wave super- 
conductors as the Holstein effecttlt 



164 





FIG. 16. a) The exchange diagram for a boson mediated 
interaction. The solid line stands for a propagating fermion. 
The wiggly line is a phonon propagator in the case of elec- 
tron- phonon interaction, and a magnon line in the case of 
a spin- fluctuation mediated interaction, b) The lowest or- 
der diagram for the fermionic self energy due to a direct four 
fermion interaction, also represented by a wiggly line. The 
flgure is taken from 149 



Consider a clean s— wave superconductor, and suppose 
that the residual interaction between fermions occurs via 
the exchange of an Einstein phonon. Assume for simplic- 
ity that the fully renormalized electron phonon coupling 
is some constant ^ep, and that the phonon propagator 
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D{q,Lo) is independent of momentum q and has a sin- 
le at a phonon frequency wq (the Holstein model) 
Phonon exchange gives rise to a fermionic self- 
energy (see Fig |l^a) 



E(ia;„ 



{2tj:Y 



G]^{iuin)D{iuJrn - i^^n) (78) 



f) 



which is a convolution of D{lo) = \/{loq — (w + iO 
with the full fermionic propagator Gk(i-ij), which in a su- 
perconductor is given by EqJlJ: 



Gk(a;) 



(79) 



As before, <^{lo) is the pairing vertex, and is the band 
dispersion of the fermions. At T = both S"(u;) and 
$"(a;) obviously vanish for w < A. This implies that the 
fermionic spectral function Ak(w) = |Gj^(w)| /tt for parti- 
cles at the Fermi surface (k = kp) has a 5— function peak 
at w = A, i.e. A is a sharp gap at zero temperature. 
The fermionic density of states in a superconductor 



N{u)) = No Im 



E(^) 



($2(^)_I]2(^))l/2 



(80) 



vanishes for lo < A and has a square-root singularity 
N(uj) (X (io — A)~^/^ for frequencies above the gap (A^o 
is the normal state density of states) . 

The onset of the imaginary part of the self-energy, 
Eq.[78|, can be easily obtained by using the spectral repre- 
sentation for fermionic and bosonic propagators in Eq. |7^ 
and re-expressing the momentum integration in terms of 
an integration over e^. At T = we obtain 



E"(w > 0) cx / duj'N{uj')D"{uj - u') 



(81) 



Since for positive frequencies, D"{lo) — {ttDo/2ljJq)5{uj - 
L0o)^ the frequency integration is elementary and yields 



Y."{uj > 0) cx N{uj-ujq). 



(82) 



We see that the single particle scattering rate is di- 
rectly proportional to the density of states shifted by the 
phonon frequency. Clearly, the imaginary part of the 
fermionic self-energy emerges only when uj exceeds the 
threshold 



rjt = A + cjo, 



(83) 



the sum of the superconducting gap and the phonon 
frequency. Right above this threshold, S"(a;) cx (w — 
ilt)"^^^- By the Kramers-Kronig relation, this non- 
analyticity causes a square root divergence of S'(w) at 
w < ilf Combining the two results, we find that near 
the threshold, S(w) = A + C/V^^t — w where A and C 
are real numbers. By the same reasoning, the pairing 



vertex $(a;) also possesses a square-root singularity at 
Vlt. Near uj = Vtt, <i>(w) = B + C/^J^lt - oj with real B. 
Since Vlt > A, we have A > B. 

The singularity in the fermionic self-energy gives rise 
to an extra dip- hump structure of the fermionic spectral 
function at k = \ip. Below fit, the spectral function 
is zero except for lo = A, where it has a 5— functional 
peak. Immediately above 17t , A{uj) (x lm{J^{uj)/ 
^^(w))) takes the form A{uj) cx (w - iltY^"^. At larger 
frequencies, A(uj) passes through a maximum, and even- 
tually vanishes. Adding a small damping introduced by 
either impurities or finite temperatures, one obtains the 
spectral function with a peak at w = A, a dip at a; « Jit, 
and a hump at a somewhat larger frequency. This be- 
havior is shown schematically in Fig. p7| . 




FIG. 17. The schematic form of the quasiparticle spectral 



function in an s— wave superconductor. Solid linq 
dashed line - at finite T. fit = A + uo (from Ref. 



T = 0, 
)• 



The singularities in S(a;) and ^{uj) affect other ob- 
servables such as the fermionic DOS, optical conductiv- 
ity, Raman respoiise, and the SIS tunneling dynamical 
conductancet£SE£2. 

For a more complex phonon propagator, which de- 
pends on both frequency and momentum, the singular- 
ities in the fermionic self-energy and other observables 
are weakeXjanji may only show up in the derivatives over 
frequencjioEEZl. Still, the opening of the new relaxational 
channel at Qt gives rise to singularities in the electronic 
properties of a phonon-mediated s— wave superconduc- 
tor. 



B. Similarities and discrepancies between d- and 
s— wave superconductors 

As we already discussed, for magnetically mediated 
d— wave superconductivity, spin fiuctuations play the role 
of phonons. Below T^, spin excitations are propagating, 
magnon-like modes with the gap A^. This A^ obviously 
plays the same role as wq for phonons, and hence we ex- 
pect that for spin-mediated pairing, the spectral function 
should display a peak-dip-hump structure as well. Fur- 
thermore, we will demonstrate below that for observables 
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such as the DOS, Raman intensity and the optical con- 
ductivity, which measure the response averaged over the 
Fermi surface, the angular dependence of the c?— wave 
gap A{9) oc cos (26*) softens the singularities, but does 
not wash them out over a finite frequency range. Indeed, 
we will find that the positions of the singularities are not 
determined by some averaged gap amplitude but by the 
maximum value of the d— wave gap, Amax = A, i.e., the 
Holstein effect is still present for a d— wave superconduc- 
tor. 

Despite many similarities, the feedback efi'ects for 
phonon-mediated s— wave superconductors, and magnet- 
ically mediated d— wave superconductors are not equiva- 
lent as we now demonstrate. The point is that for s— wave 
superconductors, the exchange process shown in Fig.[l^a 
is not the only possible source for the fermionic decay: 
there exists another process, shown in Fig, 16 d, in which 
a fermion decays into three other fermions. Th|i(S j^rocess 
is due to a residual four-fermion interaction 
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One 

can easily make sure that this second process also gives 
rise to the fermionic decay when the external uj exceeds a 
minimum energy of 3 A, necessary to pull all three inter- 
mediate particles out of the condensate of Cooper pairs. 
At the threshold, the fermionic spectral function is non- 
analytic, much like that found at A + uiq- This implies 
that in s-wave superconductors, there are two physically 
distinct singularities, at A + ujq and at 3 A, which come 
from different processes and therefore are independent 
of each other. Which of the two threshold frequencies is 
larger depends on the strength of the coupling and on the 
shape of the phonon density of states. At weak coupling, 
is exponentially larger than A, hence the 3 A thresh- 
old comes first. At strong coupling, u)q and A are compa- 
rable, but calculations within the Eliashberg formalism 
show that for real-^aaterials ( e.g. for lead or niobium), 
still 3 A < A + u)Q^2^. This result is fully consistent with 
the photoemission data for these materials 
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For magnetically mediated d-wave superconductors the 
situation is different. As we discussed in Section 2, in the 
one-band model for cuprates, which we adopt, the under- 
lying interaction is a Hubbard-type four-fermion interac- 
tion. The introduction of a spin fluctuation as an extra 
degree of freedom is just a way to account for the fact 
that there exists a particular interaction channel, where 
the effective interaction between fermions is the strongest 
due to the proximity to a magnetic instability. This im- 
plies that the spin fluctuation propagator is made out of 
particle-hole bubbles like those in Fig.[l^. Then, to the 
lowest order in the interaction, the fermionic self-energy 
is given by the diagram in Fig.p^. Higher-order terms 
convert a particle-hole bubble in Fig.^Lqb. into a wiggly 
line, and transform this diagram into the one in Figjlqa. 
Clearly then, inclusion of both diagrams would be dou- 
ble counting, i.e., there is only a single process which 
gives rise to the threshold in the fermionic self-energy. 
Note also that the fact that the diagram in Fig ^6|b is 
a part of that in Fig.p^ implies that the development 
of a singularity in the spectral function at a frequency 



different from 3A cannot be due to effects outside the 
spin-fermion model. Indeed, we will show that the model 
generates two singularities: at 3A, and at A 4- As < 3A. 
The fact that this is an internal effect, however, implies 
that As depends on A. The experimental verification 
of this dependence can then be considered as a "finger- 
print" of the spin-fluctuation mechanism. Furthermore, 
as the singularities at 3A and A -I- As are due to the 
same interaction, their relative intensity is another gauge 
of the magnetic mechanism for the pairing. We will ar- 
gue below that some experiments on cupratesp-^articu- 
larly measurements of the optical conductivityll£3, allow 
one to detect both singularities, and that their calculated 
relative intensity is consistent with the data. 

We now discuss separately the behavior of the elec- 
tronic spectral function, the density of states, SIS tun- 
neling, the Raman intensity and the optical conductiv- 
ity. To account for all features associated with d— wave 
pairing, we will keep the momentum dependence of the 
fermionic self-energy and the pairing vertex on momenta 
along the Fermi surface, although this dependence is in- 
deed weak near hot spots. For simplicity, we assume a cir- 
cular Fermi surface. In this situation, the fc-dependence 
of the self-energy and the pairing vertex reduces to the 
angular dependence, i.e. E = T,{9,uj) and <& = $(0,a;) 



C. The spectral function 

We first consider the spectral function A\f^{u}) = 
{I / ■k)\G'{^{uj)\. In the superconducting state, for quasi- 
particles near the Fermi surface 



^k(w > 0) = - Im 

TT 



E(0,c^)+ek 



{uj + T.{e,uj))-^^{e,uj) 



(84) 



By definition, A]f^{—uo) = A\^{uo). 

In a Fermi gas with d-wave pairing, Y,{9,uj) = 0, and 
<^{9,uj) = A{9) cx cos (20). The spectral function then 
has a 5— function peak &t lj — (A^(0) -I- e^)^^^- It is ob- 
vious, but essential for comparison with the strong cou- 
pling case, that the peak disperses with k and that far 
away from the Fermi surface one recovers normal state 
dispersion. 

For strong coupling we consider the spectral function, 
Ak(w), for fermions located near hot spots, 9 = 9hs- 
where the gap A{9) ( defined as a solution of Yi'iu) = 
A,9hs) = ^'{cl! — A,9hs)) is maximum. As discussed 
above, we expect the spectral function to possess a peak 
at w = A and a singularity at w = f^t = A -I- As. The 
behavior of A(io) near the singularity is robust and can 
be obtained without a precise knowledge of the frequency 
dependence of S(w) and $(aj). All we need to know is 
that near a; = A, S^(w) — ^^{uj) cx w — A. Substituting 
this form into Eq|7^ and converting to the real axis using 
the spectral representation, we obtain for uj = fit + e 
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dx 



(x(e-a;))V2 



(85) 



This integral is the same as in Eq.|7^, hence E" under- 
goes a finite jump ai uj — fit, just as the spin polarization 
operator does a.t uj — 2 A. By the Kramers-Kronig rela- 
tion, this jump gives rise to a logarithmic divergence of 
S'. The same singular behavior holds for the pairing ver- 
tex with exactly the same prefactor in front of the 
logarithm. The last result implies that — ^(lu) is 
non-singular at u) = fit- Substituting these results into 
Eq.p^, we find that the spectral function A{lj) behaves 
at w > fit as l/log^(u; — fit), i.e., almost discontinuously. 
Obviously, at a small but finite T, the spectral function 
should have a dip very near uj ^ Qt, and a hump at a 
somewhat higher frequency. 
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FIG. 18. Upper panel: the quasiparticle spectral function 
determined by solving the coupled Eliashberg equations for 
A = 1. The peak-dip-hump structure of A{uj) is clearly visible 
but not dramatic. Middle panel: real and imaginary parts of 
the fermionic self-energy (dashed and solid lines, respectively) 
. Lower panel -the frequency derivative of E"(a;). The extra 
struc ture at 3 A is clearly visible. The figure is taken from 
Ref. 
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In Fig.|T^ we present the result for A{lu) obtained from 
a solution of the-set of three coupled Eliashberg equa- 
tions at T ^ Tcllla. This solution is consistent with our 
analytical estimate. We clearly see that the fermionic 
spectral function has a peak-dip-hump structure, and the 
peak-dip distance equals Ag. We also see in FigjlJ that 
the fermionic self-energy is non-analytic at w = 3A. As 
we discussed above, this last non-analyticity originates in 
the non-analyticity of the dynamical spin susceptibility 
at w = 2A. 



Another "fingerprint" of the spin-fluctuation scatter- 
ing can be found by studying the evolution of the spec- 
tral function as one moves away from the Fermi surface. 
The argument here goes as follows: at strong coupling, 
where A > ajgf, probing the fermionic spectral function 
at frequencies progressively larger than A, one even- 
tually probes the normal state fermionic self-energy at 
UJ Wsf. Substituting the self energy Eq.^ into the 
fermionic propagator, we find that up to a; lj, the 
spectral function in the normal state does not have a 
quasiparticle peak at a; = ek- Instead, it only displays 
a broad maximum at uj — e'^/uJ- In other words, at 
ujsf < UJ < UJ, the spectral function in the normal state 
displays non-Fermi liquid behavior with no quasiparticle 
peak (see Fig. ||). The absence of a quasiparticle peak 
in the normal state implies that the sharp quasiparticle 
peak that we found at = A for momenta at the Fermi 
surface cannot simply disperse with k, as it does for non- 
interacting fermions with a d-wave gap. Specifically, the 
quasiparticle peak cannot move further in energy than 
A-|- As since at larger frequencies, spin scattering rapidly 
increases, and the fermionic spectral function should dis- 
play roughly the same non-Fermi-liquid behavior as in 
the normal state. 
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FIG. 19. a) Frequency dependence of the spectral function 
in the superconducting state for different efc. The curve at 
the bottom has a highest e^. No coherent quasiparticle peak 
occurs for energies larger than A -I- A^. Instead, the spectral 
function displays a broa d m aximum, similar to that in the 
normal state. (From Ref 149.) 



In Fig, 19 a we present plots for the spectral function 
as the momentum moves away from the Fermi surface. 
We see the behavior we just described: the quasiparticle 
peak does not move further than A -t- Ag. Instead, when 
k — kp increases, it gets pinned at A -I- As and grad- 
ually looses its spectral weight. At the same time, the 
hump disperses with k and for frequencies larger than 
A + As gradually transforms into a broad maximum at 
UJ — e\luJ- The positions of the peak and the dip versus 
k — kp are presented in Fig.n£^ 
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D. The density of states 

The quasiparticle density of states, N{u!), is the mo- 
mentum integral of the spectral function: 



(86) 



Substituting Ak(ti^) from Eq.^ and integrating over e^, 
one obtains 



N{u;) (X Im / d9 



($2(f?,co)-(t^ + S(0,w)f)i/2' 

(87) 
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FIG. 20. Density of states of a noninteracting Fermi gas 
with d-wave gap|-(solid line) and with s-wave gap (dashed 
line). (From Ref.tia). 

We first consider ]>S{ijj) in a d-wave gas, and then dis- 
cuss strong coupling effects. In a d— wave gas, E = 
and Alt [p A cos (26*). Integrating in Eq.B^ over B we 
obtaintBEB 



7V(a;) = TVo Re 



2^ 



2ir 



d6i 





2Na ( K{A/uj) 
n \{uj/A)K{uj/A) 



- A2 cos2(26i) 

for Lj > A; 
for (jj < A. 



where i4'(a;) is the elliptic integral of first kind. We see 
that N{llj) for <C a and diverges logarithmic as 
(I/tt) ln(8A/|A — for w « A. At larger frequencies, 
N{lu) gradually decreases towards the frequency indepen- 
dent, normal state value of the DOS, which we have nor- 
malized to unity. The plot of N(uj) in a d— wave BCS 
superconductor is presented in Figj20[ 

For comparison, in an s-wave superconductor, the DOS 
vanishes at w < A and diverges as N{uj) (x {uj — A)~^^^ 
at w > A. We see that a d-wave superconductor is dif- 
ferent in that (i) the DOS is finite down to the smallest 
frequencies, and (ii) the singularity at a; = A is weaker 
(logarithmic). Still, however, N{uj) is singular only at a 



frequency which equals to the largest value of the d— wave 
gap. This illustrates a point made earlier: the angular 
dependence of the d-wave gap reduces the strength of the 
singularity at uj = Ajnax{9) , but does not wash it out 
over a finite frequency range. 
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FIG. 21. (a) The behavior of the SIN tunneling conduc- 
tance (i.e., DOS) in a strongly coupled d-wave superconduc- 
tor. Main pictures - N{uj), insets - dN (uo) / dw . (a) The 
schematic behavior of the DOS for a flat gap. (b) The solution 
of the Eliashb erg- type equations for a flat gap. The shaded 
regions are the ones in which the flat gap approximation is 
incorrect as the physics is dominated by nodal quasiparticles. 
(c) The schematic behavior of N{uo) for the quadratic varia- 
tion of the gap near its maxima, (d) The expected behavior 
of the DOS in a real situation when singularities are softened 
out by finite T or impurity scattering. The position of A -I- As 
roughly correspon ds to a minimum of dN{u)/duo. The flgure 
is taken from Ref. 149. 



We now turn to strong coupling. We first demonstrate 
that the DOS possesses extra peak-dip features, associ- 
ated with the singularities in E(ti;) and ^[uj) at uj — 
which for spin-mediated pairing is 



(89) 



An analytical approach proceeds as followalZ3. Consider 
first a case when the gap is totally fiat near a hot spot. 
At w = 17t, both S(a') and ^{iS) diverge logarithmically. 
On substituting these forms into Eq. we find that 
N{ijj) has a logarithmic singularity: 



iVsing(w) OC log 



1 



1/2 



(90) 



This singularity gives rise to a strong divergence of 
dN(uj)/du) at w = fit, a behavior schematically shown in 
Fig. (|2l|a) . In part (b) of this figure we present the result 
for AYcij) obtained by the solution of the Eliashberg-type 
Eqs.pqj72. A small but finite temperature was used to 
smear out divergences. We recall that the Eliashberg set 
does not include the angular dependence of the gap near 
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hot spots, and hence the numerical result for the DOS in 
Fig. pi| b should be compared with Fig. pT|a. We clearly 
see that N{uj) has a second peak at w = fit- This peak 
strongly affects the frequency derivative of N{lu) which 
become singular near Qf 

The relatively small magnitude of the singularity 
in N(uj) is a consequence of the linearization of the 
fermionic dispersion near the FermL-swrface. For an ac- 
tual eic chosen to fit ARPES data liHa, nonlinearities in 
the fermionic dispersion occur at energies comparable to 
rif This is due to the fact that hot spots are located close 
to (0, tt) and related points at which the Fermi velocity 
vanishes. As a consequence, the momentum integration 
in the spectral function should have a less pronounced 
smearing effect than found in our calculations, and the 
frequency dependence of iV(a;) should more resemble that 
of A{uj) for momenta where the gap is at maximum. 

For a momentum dependent gap, the behavior of 
fermions near hot spots is the same as when the gap 
is flat, but now fit depends on 6 as both A and As vary 
as one moves away from a hot spot. The variation of 
A is obvious, the variation of A^ is due to the fact that 
this frequency scales as A^/^. Since both A and A^ are 
maximal at a hot spot, we can model the momentum 
dependence by replacing 



fit 



(91) 



where 6 = 6 — 9hs, and a > 0. The singular pieces of 
the self-energy and the pairing vertex then behave as 
log(f2t — w — a9^)^^. Substituting these forms into Eq.|87| 
and using the fact that S(cli) — $(a;) « const at cj « fit, 
we obtain 



the angular dependence of the gap, but still persists at a 
particular frequency related to the maximum value of the 
gap. This point is essential as it enables us to read off the 
maximum gap value directly from the experimental data 
without any " deconvolution" of momentum averages. 

For real materials, in which singularities are removed 
by e.g., impurity scattering, N{uj) likely has a dip at 
Lu > rit and a hump at a larger frequency. This is 
schematically shown in Fig. ([2l|)d. The location of Qt is 
best described as a point where the frequency derivative 
of the DOS passes through a minimum. The singular- 
ity in dN[ijj)/dLO at f2t, and the dip-hump structure of 
N{bj) at cj > fit are additional "fingerprints" of the spin- 
fluctuation mechanism in the single particle response. 



E. SIS tunneling 



Measurements of the dynamical conductance dl/dV 
through a superconductor - insulator - superconductor 
(SIS) junction offer another tool to search for the finger- 
prints of the spin-fluctuation mechanism. The conduc- 
tance through this junction is the dqiisiative over voltage 
of the convolution of the two D0sll2il: dl/dV oc S{lo) 
where 



dnN{uj-n) dnN{n). 



(94) 



N.ingiio) (xRe I d0 log(17t 



a6l^ 



-1/2 
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A straightforward analysis of the integral shows that now 
N{uj) displays a one-sided non-analyticity at a; = fit: 



iVsi„g(w) = -Be{Lu~ nt) 



log(w -fit)] 



1/2 



(93) 



where B > 0, and 8(x) = 1 for x > 0, and 8(x) = for 
X < 0. This non-analyticity gives rise to a cusp in N(uj) 
right above f2t, and one-sided square-root divergence of 
the frequency derivative of the DOS. This behavior is 
shown schematically in Fig. Comparing this be- 

havior with that shown in Fig. pT| a for a flat gap, we 
observe that the angular dependence of the gap predom- 
inantly affects the form of N{u!) at w < fit- At these 
frequencies, the angular variation of the gap completely 
eliminates the singularity in N{uj). At the same time, 
above fit , the angular dependence of the gap softens the 
singularity, but, still, the DOS sharply drops above fit in 
such a way that the derivative of the DOS diverges on ap- 
proaching fit from above. We see again that in a d— wave 
superconductor, the singularity in the DOS is softened by 



The DOS in a d— wave gas is given in Eq. 88. Substitut 



ing this form into Eq. p4 and integrating over frequency 
we obtain the resulLpresented in Fig|2^. At small w, <S'(i^) 
is quadratic in aJl£3. This is an obvious consequence of 
the fact that the DOS is linear in us. At w = 2A, S{oj) 
undergoes a finite jump. This jump is related to the 
fact that near 2A, the integral over the two DOS in- 
cludes the region around = A where both N{fl) and 
N{u! — fl) are logarithmically singular, and 9oiV(r2) di- 
verges as 1/(17 — A). The singular contribution to S{lu) 
from this region can be evaluated analytically and yields 



1 



-P 



dx log|a;| 



1 



sign(w — 2A) 



(95) 



Observe that the amplitude of the jump in the SIS con- 
ductance is a universal number which does not depend 
on the value of A. At larger frequencies, S{uj) contin- 
uously goes down and eventually approaches a value of 
S{uj oo) — 1. 



29 



S,(0, a) -)*^^ 




FIG. 22. (a) The schematic behavior of the SIS tunnehng 
conductance, S{uj), in a strongly coupled d-wave supercon- 
ductor. Main pictures - ^(tj), insets - dS{u}) / duj. (a) The 
schematic behavior of S{uj) for a flat gap. (b) The solution 
of the Eliashb erg- type equations for a flat gap using the DOS 
from the previous subsection. The shaded regions are the ones 
in which the flat gap approximation is incorrect as the physics 
is dominated by nodal quasiparticles. (c) The schematic be- 
havior of S{(jj) for a quadratic variation of the gap near its 
maxima, (d) The expected behavior of the SIS conductance 
in a real situation when singularities are softened out by finite 
T or by impurity scattering. 2A -I- As roughly corresponds to 
the maximum of dS{u)) / duo. (From Ref. 14t). 



In the case of strong coupling one finds again that the 
quadratic behavior at low frequencies and the disconti- 
nuity at 2A survive at arbitrary coupling. Indeed, the 
quadratic behavior at low uj is just a consequence of the 
linearity of N{u)) at low frequencies. Therefore, just as 
we did for the density of states we concentrate on behav- 
ior above 2A that is sensitive to strong coupling effects. 

Consider first how the singularity in S(tj^) at f2t affects 
S{uj). From a physical perspective, we would expect a 
singularity in S{ijj) at cj = A -|- Jit = 2A + Aj. Indeed, 
to get a nonzero SIS conductance, one has to first break 
a Cooper pair, which costs an energy of 2A. After a pair 
is broken, one of the electrons becomes a quasiparticle 
in a superconductor and takes the energy A, while the 
other tunnels. If the tunneling voltage equals A -I- fit, the 
electron which tunnels through a barrier has energy Sit, 
and can emit a spin excitation and fall to the bottom of 
the band (see Fig. |2^). This behavior is responsible for 
the drop in dl/dV and is schematically shown in Fig. |23. 




FIG. 23. The schematic diagram for the dip features in 
SIN and SIS tunneling conductances (figures a and b, re- 
spectively) . For SIN tunneling, which measures the fermionic 
DOS, the electron which tunnels from a normal metal can 
emit a propagating magnon if the voltage eV = A -I- . Af- 
ter emitting a magnon, the electron falls to the bottom of 
the band. This leads to a sharp reduction of the current and 
produces a drop in dl/dV. For SIS tunneling, the physics 
is similar, but one first has to br eak an electron pair, which 
costs energy 2 A (taken from Ref. 172). 



Consider this effect in more detailEZ3Ef3. We first note 
that cj = A -I- Sit is special for Eq. |9^ because both 
dN{n) /dil and N{u! — ft) diverge at the same energy, 
i7 = ilt. Substituting the general forms of N{lu) near 
UJ = fit and w = A, we obtain after simple manipula- 
tions that for a flat gap, S'(w) has a one-sided divergence 
at w = fit + A = 2A + As. 



S'sing(e) oc 



e(-e) 
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where e = co — {fit + This obviously causes a di- 
vergence of the frequency derivative of <S'(a-') (i.e., of 
d^I/dV^). This behavior is schematically shown in Fig. 

In Fig. nib we present the results for S{oj) ob- 
tained by integrating theoretical N{u!) from Fig.plp. We 
clearly see that S{ui) and its frequency derivative are sin- 
gular at UJ = 2 A -|- As , in agreement with the analytical 
prediction. 

For a quadratic variation of the gap near the maxima, 
calculations similar to those for the SIN tunneling yield 
the result that S{uj) is continuous through 2 A + As, but 
its frequency derivative diverges as 



dSjuj) 
duj 



oc P 



dx 



e(-e) 



(a;|loga;|)i/2(a;„e) |elog |e||i/2 ' 

(97) 



The singularity in the derivative implies that near e = 

1/2 



^(e) = 5(0) - C e(-e) 



log(-e)| 



(98) 



where C > 0. This behavior is schematically presented 
in Fig. 22 i. We again see that the angular dependence 
of the gap softens the strength of the singularity, but 
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the singularity remains confined to a single frequency 
u = 2A + As. 

In real materials, the singularity in S(lu) is softened 
and transforms into a dip slightly below 2 A + A^, and 
a hump at a frequency larger than 2 A + A^. The fre- 
quency 2A + As roughly corresponds to a maximum of 
the frequency derivative of the SIS conductance. 



F. Optical conductivity and Raman response 

Further phenomena sensitive to fit are the optical con- 
ductivity, cr(w) and the Raman response, R{lo). Both are 
proportional to the fully renormalized particle-hole polar- 
ization bubble, but with different signs attributed to the 
bubble composed of anomalous propagators. Specifically, 
after integrating in the particle-hole bubble over £k, one 
obtains 



i?(w) = im / dui'dev^ {e)u,{e,Lu,Lu') 



— Re 



(99) 



where V{0) is a Rpaan vertex which depends on the 
scattering geometryllZ3, and 



?,w,u;') 



D+D^{D+ + D^ 
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Here a = — 1 for 11^, and a = 1 for H^r- Also, S± = 
E(a;±) and $± — ^{w±), where u!± — u' ± w/2. We 
also introduced D± = ($| - S|)i/2. Note that S and $ 
depend on uj and 0. 

In a superconducting gas, the optical conductivity van- 
ishes identically for any nonzero frequency due to the 
absence of a physical scattering between quasiparticles 
in a gas. The presence of a superconducting condensate, 
however, gives rise to a S functional term in cr at u; = 0: 
(t{uj) = 7r6{uj) J d9dLLj'Ilcr{9,Q,uj'). JJJnifi behavior is typ- 
ical for any BCS superconductorEZEI. The behavior of 
a{Lu) for a d-wave gas with additional impurities, caus- 
ing inelastic scattering, is more complex and has been 
discussed in Ref. 18C. 



The form of the Raman intensity depends on the scat- 
tering geometry. For the scattering in the Big channel, 
the Raman vertex has the same angular dependence as 
the d-wave gap, i.e., ¥{9) cx cos {29ytIMB. Straightfor- 
ward computations then show that at low frequencies, 
R{lu) cx cj'^EIj. For a constant V{9), we would have 

R{lj) (X LU. 

Near to = 2A, the Big Raman intensity is singular. For 
this frequency, both and _D_ vanish at w' = and 
9 = 0. This causes the integral for R{lu) to be divergent. 
The singular contribution to R{uj) can be obtained ana- 
lytically by expanding in the integrand to leading order 
in and in 9. Using tbe-^pectral representation, we then 
obtain, for lu = 2A + il3 



R{uj) = f dVL [ d9 \ 

Jo J V^ + a9'^VS -n + a9'^ 
1 



^n + a9^ + \/S-n + a9^) 



(101) 



where, as before, 9 = 9 — 9hs For a fiat band (a = 0), 
R{uj) cx i?e[((jj — 2A)^^/^]. For a 7^ 0, i.e., for a quadratic 
variation of the gap near its maximum, the 2d integration 
in Eq.lOl is elementary, and yields i?(w) cx log \u) — 2A|. 
At larger frequencies i?(w ) gradually decreases. 

The behavior of Riyi) in a d-wave gas is shown in 
Fig.^. Observe that due to the interplay of numeri- 
cal factors, the logarithmic singularity shows up only in 
the near vicinity of 2A, while at somewhat larger w, the 
angular dependence of the gap becomes irrelevant, and 
_R(aAJpehaves as (w — 2A)^/^, i.e. the same as for a flat 
gaptlj. 
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FIG. 24. The behavior of the Raman response in a BCS su- 
perconductor with a flat gap (dashed line), and for a quadratic 
variation of the gap near its maximum (solid line). 

We now consider strong coupling effects. A nonzero 
fermionic self-energy mostly affects the optical conduc- 
tivity for the simple reason that it becomes finite in the 
presence of spin scattering which can relax fermionic mo- 
menta. For a momentum-independent gap, a finite con- 
ductivity emerges above a sharp threshold. This thresh- 
old stems from the fact that at least one of the two 
fermions in the conductivity bubble should have a fi- 
nite E", i.e., its energy should be larger than fit. An- 
other fermion should be able to propagate, i.e., its energy 
should be larger than A. The combination of the two re- 
quirements yields the threshold for aito > 0) at 2A-l-ilt, 
i.e., at the same frequency where the SIS tunneling con- 
ductance is singular. One can easily demonstrate that for 
a fiat gap, the conductivity behaves above the threshold 
as e^/^/ log^ e, where e = uj - {A + = lu - (2A + A^). 
This singularity obviously causes a divergence of the first 
derivative of the conductivity at e = -1-0. 

In Fig.^ we show the result for the conductivity ob- 
tained by solving the set of coupled Eliashberg-type equa- 
tions, Eqsj6^[72 109,181. We see the expected singularity 



at 2 A + As. The insert shows the behavior of the inverse 
conductivity 1/(7{lu) Observe that ^1/^7(0;) is linear in lo 
over a rather wide frequency rangeE 
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FIG. 25. The real part of the optical conductivity (Ji(tj) at 
the lowest T obtained using the self energy and the pairing 
vertex from the solution of the Eliashberg equations for A = 1. 
The onset of the optical response is = 2 A + A^. The con- 
tributions from nodal regions (not included in calculations) 
yield a nonzero conductivity at all uj. Inset - the behavior of 
the inverse conductivity vs frequency. (From Ref. 



109) 



For a true d— wave gap, the conductivity is finite for 
all frequencies simply because the angular integration in 
Eq. ^ involves the region near the nodes, where S" is 
nonzero down to the lowest frequencies. Still, the conduc- 
tivity is singular at f2t + A as we now demonstrate. In- 
deed, as we already discussed, at deviations from 9 = dhs, 
where the gap is at maximum, both A and As decrease, 
hence ftt{d) = fit — o-9^, where 9 = 9 — 9hs and a > 0. 
The singular pieces in S(w) and ^{lo) then behave as 
I log(rit — w — a9'^)\. Substituting these forms into the 
particle-hole bubble and integrating over 9, we find that 
the conductivity and its first derivative are continuous at 
UJ = 2 A -I- As, but the second derivative of the conduc- 
tivity diverges as (Pa/dw'^ oc l/(|e| log^ e). We see that 
the singularity is weakened by the angular dependence of 
the gap, but is still located exactly at rit + A = 2A-|-As. 

The same reasoning can be applied to a region near 
4 A. The singularity at 4 A is also weakened by the an- 
gular dependence of the gap, but is not shifted and still 
should show up in the second derivative of the conduc- 
tivity. 

For the Raman intensity, strong coupling effects are 
less relevant. First, one can prove along the same lines 
as in previous subsections that the cubic behavior at low 
frequencies for Big scattering (and the linear behavior 
for angular independent vertices), and the logarithmic 
singularity at 2A, are general properties of a d— wave 
superconductor, which survive for all couplings. Thus, 
analogous to the density of states and the SIS-tunneling 
spectrum, the Raman response below 2A is not sensitive 
to strong coupling effects. Second, near loq + A, sin- 
gular contributions which come from and ^-i,^- 
terms in 11^ in Eq. |9^ cancel each other. As a result, for 
a flat gap, only the second derivative of R{uj) diverges 
at A -|- rif For a quadratic variation of a gap near its 
maximum, the singularity is even weaker and shows up 
only in the third derivative of R{lo). Obviously, this is a 
very weak effect, and its experimental determination is 
difficult. 



We now argue that measurements of the optical con- 
ductivity allow one not only to verify the magnetic sce- 
nario, but also to determine both As and A indepen- 
dently in the same experiment. In the magnetic scenario, 
the fermionic self-energy is singular at two frequencies: 
at fit = A -|- As, which is the onset frequency for spin- 
fluctuation scattering near hot spots, and at w = 3 A, 
where fermionic damping near hot spots first emerges 
due to a direct four-fermion interaction. Since in the 
spin-fluctuation mechanism, both singularities are due to 
the same underlying interaction, their relative intensity 
can be obtained within the model. 

In general, the singularity at 3A is much weaker at 
strong coupling, and can be detected only in the analysis 
of the derivatives of the fermionic self-energy. We recall 
that the singularity in S(ti') at fit gives rise to singularity 
in the conductivity at A -f fit, while the 3 A singularity 
in S(a;) obviously causes a singularity in conductivity at 
Lo = 4A. In addition, we expect a singularity in cr(w) 
at 2f2t, as at this frequency both fermions in the bubble 
have a singular S(rit)- 
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FIG. 26. The calculated frequency dependence of 
W(ijj) = ^[cjRe[l/cr((.j)]] at r ^ 0. This quantity is a sen- 
sitive measure of the fine structure in the optical response. 
The locations of the extrema are: (1) 2A -I- A^, (2) 4A, (3) 
2 A -|- 2 As . Observe that the maximum shifts to a lower tem- 
perature, but the minimum remains at 2A -f- As . (From Ref. 

looj) 



For superconductors with pairing due to electron- 
phonon interaction the fine structure of the optical con- 
ductivity has been analyzed by studying the second 
derivative of conductivity via Wiui) = ^p^(u;Re(T~^(w)) 
which is proportional to c?{uj)F(uj) where aiuS) is an ef- 
fectiye-electron-phonon coupling, and F(ui) is a phonon 
DOSllij. In Figj2g we presenL-the result of the strong 



coupling calculations of 1/^(^^)1123. There is a sharp max- 
imum in W{uj) near 2A + As, which is followed by a deep 
minimum. This form is consistent with our analytical ob- 
servation that for a flat gap (which we used in our numer- 
ical analysis), the first derivative of conductivity diverges 
at = 2 A -|- As. At a finite T (a necessary attribute of a 
numerical solution), the singularity is smoothed, and the 
divergence is transformed into a maximum. Accordingly, 
the second derivative of the conductivity should have a 



32 



maximum and a minimum near 2 A + Ag. The numeri- 
cal analysis shows that the maximum shifts to lower fre- 
quencies with increasing T, but the minimum moves very 
little from 2A -|- As, and is therefore a good measure of 
a magnetic "fingerprint" . 

Second, we note from Fig.|2^ that in addition to the 
maximum and the minimum near 2 A -I- Ag, W{uj) has 
extra extrema at 4 A and 2^lt = 2 A -I- 2 As. These are 
precisely the extra features that we expect: they are a 
primary effect due to a singularity in T,(uj) ai uj — 3A and 
a secondary effect due to a singularity in at = lit- 
The experimental discovery of these features will be a fur- 
ther argument in favor of spin-mediated pairing and the 
applicability of the spin-fermion model. 



VI. COMPARISON WITH THE EXPERIMENTS 
ON CUPRATES 

In this section we compare the theoretical results 
for the spin-fermion model of the nearly antiferromag- 
netic Fermi liquid with the experimental data for op- 
timally doped members of the Bi2Sr2CaCu208 and 
YBa2Cu307_j^ families of cuprate superconductors. We 
make the assumption that at this doping level, ab- 
sent the influence of a superstructure and imper- 
fections that make the NMR lines in Bi2Sr2CaCu2 08 
broad and difficult to interpret, the normal state be- 
havior of Bi2Sr2CaCu208 will resemble closely that of 
YBa2Cu307_j^. This enables us to take the two input 
parameters of the model from fits to NMR in the latter 
material. We then can compare theory and experiment 
in the normal state and as T ^ in the superconducting 
state. Finally, we discuss the general phase diagram of 
the cuprates and the pseudogap physics of these materi- 
als. 



A. Parameters of the model 

The two input parameters of the theory are the cou- 
pling constant A and the overall energy scale Cu = AX^ujsf. 
Alternatively, we can re-express A as A = 3vf£^^^ / (IGlUsi) 
and use vfS.~^ and Wsf as inputs. The values of Wsf 
and ^ can be extracted from the NMR measurements 
of the longitudinal and transverse spin-lattice relaxation 
rates, and from neutron scattering data, which measure 



Lo/usiy). We will 
'-optimally doped 
i94| yields a mod- 



primarily rely on NMR data for ru 
YBa2Cu306+5. The NMR analysisB^ 
erately temperature dependent ujsf and ^ which take the 
values LOgf ^ 15 — 20meV and ^ ~ 2a in the vicinity of Tcr, 
which for slightly overdoped materials will be close to Tc- 
The neutron data from inelastic scattering (INS) exper- 
iments on the normal state are more difhcult to analyze 
because of the background which increases the measured 
width of the neutron peak and because of the possible 



influence of weak intrinsic inhompgeneities on a global 
probe such as INS. The data showtlj that the dynamical 
structure factor in the normal state is indeed peaked at 
q = Q = (7r/a,7r/a), and that the width of the peak 
increases with frequency and at lu ^ 50meV reaches 1.5 
of its value at cj = 0. A straightforward fit to the the- 
ory yields LUgf ~ 35 — 40meV and a weakly temperature 
dependent ^ ^ a which are, as expected, larger than the 
ujsi and smaller than the ^ values extracted from NMR. 
We will be using cjgf ^ 20meV and ^ = 2a for further 
estimates. 



The value of the Fermi velocity can be obtained 
from the photoemission data on Bi2Sr2CaCu208 at 
high frequencies, where the self-energy corrections to 
the fermionic dispersion become relatively minor. We 
note that because of problems related to the sur- 
face reconstruction in YBa2Cu306+5 the vast ma- 
jority of high quality angular resolved photoemission 
spectroscopy (ARPES) experiments are performed on 
Bi2Sr2CaCu208, the material where there are much 
less reliable NMR experiments in part because of 
superstructure induced line broadening. The three 
groups that report MDC (momentum distribution curve) 
data for, B^pSi|EC aCu208 and momenta along the zone 
'^^^ all agree that the value of the bare 
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diagonal 

Fermi velocity along the diagonal (determined at higher 
energies where mass renormalization is assumed to be 
small) is rather high: 2.5 — 3eVA, or 0.7 — O.SeVa where 
a ~ 3.8A is the Cu — Cu distance. We can use the 
t — t' tight binding model for the electronic disper- 
sion to relate this velocity with that at hot spots. Us- 
ing the experimental facts that the Fermi surface is lo- 
cated at k « (0.47r/a, 0.47r/a) for momenta along the 
zone diagonal and at k w ("""/a, 0.27r/a) for along the 
zone boundary, we find t ~ 0.2 — 0.25eV, t' w —0.35t 
and /i ~ — l.li. These numbers agree with those used 
in numerical studiesli£ZI. The hot spots are located at 
khs = (O.IGtt, 0.847r) and symmetry related points, and 
the velocity at a hot spot is then approximately half of 
that along zone diagonal. This yields vp ~ 0.35 — 0.4eVa. 



Combining the results for vp, C and Wgf, we obtain 
A 1.5 - 2. This in turn yields iIi - 0.2 - 0.3eV. As 
an independent check of the internal consistency of these 
estimates, we compare theoretical and experimental val- 
ues of the resonance spin frequency As- As we said at 
the end of Section. 4, A^ ~ 0.2(D for A = 2. Substituting 
the value of a), we obtain A^ close to the experimental 
value of 40meV. A smaller ujsf = 15meV would require 
a slightly larger A, but variations of this magnitude are 
certainly beyond the quantitative accuracy of our theory. 
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FIG. 27. a. The theoretical result for the effective velocity 
of the quasiparticles Vp = vf/(1 + dT,' [cj) / duj) . For defi- 
niteness we used ujsf ~ 20meV , X — 1.7 and bare velocity 
vf = 3eVA along the diagonal, b. Experimental result £e« 
the effective velocity, extracted from the MDC dispersiont^j 
along the zone diagonal. Observe the bump in the frequency 
dependence of the velocity at 70 — 80meV in the data and at 
about 3 — 4:UJsf in the theory. 

Away from hot spots, the effective coupling decreases 
as A(k) = A/(l + where Sk is the momentum 

deviation from a hot spot along the Fermi surface. The 
largest Sk^ is for k vectors along the zone diagonals. At 
optirasO, doping, ARPES data yield Sk"'"'' ~ 0.27r/a « 
0.6/all^. We see that A is reduced by at most 1.7 as one 
moves from hot spots to the zone diagonal. A predic- 
tion of the model is that ujsf{k) increases at deviations 
from hot spots. This increase, however, should be at 
least partly compensated by the fact that ajgf oc sin^o, 
where 0o is the angle between Fermi velocities at k and 
k + Q, with 00 ~ -I in the vicinity of hot spots. (/)o 
tends to TT as k approaches the zone diagonal, and this 
reduces ujsf. In view of this competing effect which we 
cannot fully control, we believe that the effective Wsf(k) 
can best be obtained from the fit to the photoemission 
data, particularly from the MDC measurements of the 
electronic dispersion uj + = eic- In Fig ^ we com- 

pare our (1 + dY,' [oj) / doj) with the measured variation 
of the effective velocity fWw) of the electronic disper- 
sion along zone diag onalE3. We see that the theoretical 
dispersion has a bump at a; ~ Swsf (kdiag). The exper- 
imental curve;s_l{jialijpC|JLiite similar and show a bump at 
1 This yields a;sf (kdiag) ~ 25meV, 
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a value only slightly larger than that near hot spots. 

Note in passing that although 5k(^ does not vary much 
when k moves along the Fermi surface, the fact that 
the Fermi velocity is fairly large implies that along the 
zone diagonal, ekp-i-Q is roughly \/2v'pQ.2t: / a ~ 0.8eT^, 
i.e., it is comparable to the bandwidth. This implies 
that the Fcrmi-surface is very different from the near- 
perfect square that one would obtain for only nearest 
neighbor hopping. Furthermore, the fact that the Fermi 
velocity is large implies the physics at energies up to 
few hundred meV is confined to the near vicinity of 
the Fermi surface, when one can safely expand to 
linear order in k— kp. Finally, van-Hovp ^insularities 



(which we neglected) do play some roleE£3E2il but as 



Hove singularity softens due to fermionic incoherence and 
should not substantially affect the physics. The value of 
£(o.ir/a) might however be affected by an additional bi- 
layer splitting which moves one of the bands closer to 
k = (0,7r/a). 

Finally, in the analysis of the spin-fermion model we 
have neglected the temperature dependence of the cor- 
relation length, and thus of ujsi- Fits to NMR ex- 
periments on the near optimally doped member of the 
YBa2Cu306+5 family show that at Tcr — both cjsf and 
^ display mean field behavior with ^"^ ~ ^"^ (1 + T/Tq) 
and Wsf^^ — 70meV almost independent on T. 

From a theoretical perspective, the leading tempera- 
ture dependence of ^ arises from an interaction between 
spin-fluctuations and near the critical point in two dimen- 
sions has the form ^-^(T) = ^-^(T = 0) + Br4riogT 
where B = 0(1), and r4 is the effectyp-fiooir boson inter- 
action that is made out of fcrmionsEZrEj'tlj. The TlogT 
factor is the universal contribution from the bosonic loop, 
confined to momenta near Q. The four-boson interac- 
tion has two contributions: one comes from low-energy 
fermions and is universal; the other comes from high- 
energy fermions and-depends on the fermionic bandwidth 
W . One can showlllil that the temperature correction 
to ^ involves only the high-energy part of the interac- 
tion (this is what we labeled as r4). The magnitude 
of r4 can be easily estimated to be g'^ /W^. Accord- 
ingly, the temperature correction to ^ scales as '^{g/WY- 
As we have remarked, the theory is universal as long 
as 5 ^ W . In this limit, the temperature dependence 
of ^ is obviously small and can be neglected. This is 
what we will do. Notice however that in the oppo- 
site limit, when g 3> W ^ the full four-boson interaction 
differs from the lowest order term in g and might be 
estimated within an RPA-type summation. Estimates 
show that in this limit, the full r4 is fully determined 
within the low-energy sector and scales as 0(1/ J) where 
J ~ W^/g is the magnetic exchange integral. This in 
turn yields a much stronger temperature dependence of ^: 
(-^{T)-^-^{T = 0) - (T/J) logT. This resuh is similar 
to that obtained using a panlmpar cr— model approach to 



near antiferromagnetismc^ — . The agreement becomes 
obvious in the limit of a large spin-fermion interaction 
(which, we recall, is the Hubbard U if we derive the spin- 
fermion model within the RPA); double occupancy is en- 
ergetically unfavorable and the spin susceptibility obeys 
the constraint J (Pqduj xi'ij^^) c>c 1 — a;. This is equiva- 
lent to imposing a constraintj-on jthe length of the spin 
field in the nonlinear cr-modelE23 cBi . 



-{0,ir/a) 



0.34t ~ 85meV » Wsf , we expect that the van- 



B. The normal state 

In this section we compare the experimental and theo- 
retical forms of the fermionic spectral function and opti- 
cal conductivity in the normal state. No free parameters 
remain, since those which are needed to specify the model 
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completely have been taken from NMR and ARPES ex- 
periments. The discussion will follow Ref. fill. 



1 . The spectral funetion: 

The quasiparticle spectral function at various mo- 
menta is measured in angle resolved photoemission ex- 
periments. In a sudden approximation (an electron, 
hit by light, leaves the crystal without further inter- 
actions with other electrons and without paying at- 
tention to selection rules for the optical transition to 
its final state), the photoemission intensity is given by 
^k(w) = A-k^{uj)np{u!) where np is the Fermi function 
and A\^{il) ~ (I/tt) |ImG(k, is the spectral function. 
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FIG. 28. A comparison between the theoretical results of 
the spin fermion model and the photoemission MDC data. 
For the Lorentzian line-shape of the MDC signal, observed in 
experiments, the MDC hwhm equals to T," /vp- Upper panel 
- the results for the MDC hwhm vs freq uenc y at a given T. 
The experimental results are taken from 155 . Lower panel - 
the MDC fwhm vs temperature at f2 = 0. The experimen- 
tal results (righ t figure and the points oni-the left figure) are 
taken from 196, The figure is taken fromEil. 



We first use our form of the fermionic self-energy to fit 
MDC data which measure the width of the photoemission 
peak as a function of k at a given frequency. In Fig. 
we compare the theoretical results for Afc — S"(k, n)/jitp. 
with the measured Afc versus frequency at T ~ lOOK 
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We used A = 1.7 



and versus temperature at fl 
and Wgf = 20ineV. The slope of Afc is chiefly controlled 
by A. We obtain rather good agreement with the data, 
both for the frequency and temperature dependence of 
the self-energy. On the other hand, the magnitude of our 
E" is smaller than that found in the experimental data. 
To account for the values of Afc, we had to add a con- 
stant of about 70meV to S". The origin of this constant 



is unclear and explaining it is presently a challenge to 
the theorji^t may be the effect of elastic scattering by 
impuritiesE23, but the large value of this constant implies 
that it is more likely the contribution from scattering 
channels that we ignored. It is essential, however, that 
the functional dependence of Y."{il,T) can be captured 
in the spin-fluctuation approach. 

In Fig ^ we present the results for the half width at 
half maximum of the EDC (energy distribution curve) 
which measures fermionic ~ A]s_(Q)np(Q) as a 

function of frequency at a given k. For a Lorentzian line- 
shape, the EDC hwhm is given hy + 
The data are taken at T = llSKliM We see that the 
theoretical slope agrees reasonably well with the exper- 
imental one. The visible discrepancy is most likely asso- 
ciated with the fact that the experimental line-shape is 
not a perfect Lorentzian. 




FIG. 29. A comparison of the theoretical result for 
with the experimental hwhm of the EDC 
dispersion from 153 . The figure is taken from Bll 



S. The optical conductivity 



In Fig.^pSsie compare the theoretical results for the 
conductivityllHj with the experimental data for ai and 
a2 at different temperatureai^J. The theoretical results 
are obtained using the same A = 1.7 and LOgi = 20meV 
as in the fit to the photoemission data along zone di- 
agonals. Changing A affects the ratio (72 /ci at high 
frequencies, but doesjiot change the functional forms 
of the conductivitieal2j. The value of the plasma fre- 
quency was adjusted to match the d.c. conductivity and 
S" extracted from the MDC photoemission data using 
vp ^ 3eVA. This adjustment yields Wpi ~ 20000cm~^. 
This value is somewhat larger than Wpi ~ 16000cm~^ 
obtained experimen tally by integrating ai up to about 
2 — 2.5eVc3i23'E2il, however it agrees with the theoret- 
ical result 
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that the sum rule for cri(w) is exhausted 
only at extremely high frequencies of about 50a;, that are 
much larger than 2eV. 
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FIG. 30. The theoretical and experimental results for the 
real and ima ginary parts of optical conductivity. The data 
are from |l90[ The figure is taken from [HI} 

We see that theoretical calculations of Ci (oj) and (72 (oj) 
capture the essential features of the measured forms of 
the conductivities. In particular, the curves of ai at dif- 
ferent temperatures cross such that at the lowest frequen- 
cies, the conductivity decreases with T while at larger 
frequenc ies it increases with T, a behavior seen in Ref. 
195| , 191 . The imaginary part of conductivity decreases 
with T at any frequency, and the peak in cr2(w) in- 
creases iji raagnitude and shifts to lower T with decreas- 



ing 



7^(190 



191 



18£ 



At frequencies above 1500cm~^ both <Ji 
and (72 depend weakly on T and are comparable in am- 
plitude. 
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The theoretical and experimental results for 



n'j(w) = Airaiuj/ujli (from The data are from 19' 



To make the comparison more quantitative, in Fig ^ 
we present experimental and theoretical results for the 
imaginary part of the full particle-hole polarization bub- 
ble n^(w) = inaiuj/w^i. Theoretically, at T = 0, 
saturates at a value of about 0.2 independently of A and 
ppendent of frequency over a very wide 
The experimental data also clearly 



remains almost ind 



frequency range 
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show a near saturation of XI" at a value close to 0.2. 
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FIG. 32. The theoretical and experimental res ults for 
1/r* = uRea / Ima (fronfl). The data are from 190. 



The agreement between theory and experiment is, 
however, not a perfect one. In Fig. ^ we show theo- 
retical and experimental results for l/r* = U!ai/a2- The 
advantage of comparing l/r* is that this quantity does 
not depend on the unknown plasma frequency. We see 
that while both experimental and theoretical curves are 
linear in frequency, the slopes are off roughly by a factor 
of 3. This discrepancy is possibly related to the fact that 
in the spin-fermion model, W^{uj) at high enough frequen- 
cies is roughly 3 times larger than nj^(ci;)E£3, and hence 
f2/ci ^ 3, whereas experimentally cti and cr2 are com- 
parable in magnitude at high frequencies. The discrep- 
ancy in 1 /r* indicates that either the averaging over the 
Fermi surface, vertex corrections inside a particle-hole 
bubble, or RPA-type corrections to the conductivitj-T^ 
play some role. Still, Figs. ^ and |3l| indicate that the 
general trends of the behavior of the conductivities near 
optimal doping are reasonably well captured within the 
spin-fluctuation approach. 



C. The superconducting state 

In this section, we apply our results from Section 5 to 
cuprates and examine to what extent the "fingerprints" 
of spin-fluctuation pairing have been detected in experi- 
ments on optimally doped high materials. 



1. The spin susceptibility 

The major prediction of the spin fermion model for 
the spin susceptibility in the superconducting state is the 
emergence of the resonance peak in x"(Q) w) at a; = Ag if 
As < 2A. The magnitude A^ is fully determined within 
the theory and is chiefly set by the magnitude of the su- 
perconducting gap as well as the energy scale of magnetic 
fluctuations in the normal state, Wgf. For small doping 
concentration A^. cx must decrease as one approaches 
the antiferromagnetic state. The resonance mode is con- 
fined to a small region in momentum space (where it is of 
high intensity). For momenta away from Q and its close 
vicinity, magnetic excitations couple to gapless, nodal 
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quasiparticles and become overdamped, eliminating the 
resonance mode. 

In Fig.^ we show representative experimental data 
for x"(Q,cj) showing the resonance peak. at a; « 41 
meV for optimally doped YBa2Cu306.9tIll. As noted 



earlier, the position of the peak is consistent with the 
prediction of the spin fermion model. Similar behavior is 
found in Bi2Sr2CaCu208L£3; here the peak is at 43 meV. 
With, nTj^PYdoping, the measured resonance energy goes 

rderdoped YBa2Cu306.6, it 
The existence of the peak 



dowdiSiluHj. In strongly 
is approximately 25 meV 
and the downturn with underdoping agree with the pre- 
dictions of the spin-fluctuation theory. Further, the mea- 
sured amplitude of x"(Qi^) decreases, above the peak, 
but increases again for 60 — 80 meVll£3'EZ3. This might 



indeed be a 2 A effect, which appears naturally within the 
model. 



up at rct£2HlI3. All these features 
level of a weak coupling approac"^^^ 
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ly present on the 
L-Jiave been ex- 
plained within the spin-fermion modelH2l't£3. The broad.- 
ening of the peak was recently studied in detail in Ref .11^ . 
The explanation of these effects, however, requires care- 
ful analysis of the details of the electronic structure and 
is beyond the scope of this Chapter. 
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FIG. 33. Inelastic neutron scattering intensity for momen- 
tum Q = (7r,7r ) as function of frequency for YBa2Cu306.5. 
Data from Ref. 161, 



The full analysis of the resonance peak requires more 
care as (i) the peak is only observed in two-layer ma- 
terials, and only in the odd channel, (ii) the momen- 
tum dispersiaa_.of the peak is more complex than that 
for. magHjOns EZj, (iii) the peak broadens with underdop- 
ingl£2l't£3, and (iv) in underdoped materials, the peak 
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FIG. 34. ARPES spectrum for near optimally doped 
Bi22 12 for momenta close to the hot spots. Data from Ref. 
147l 



emerges at the onset of the pseudogap and only sharpens 



2. The spectral function: 

The predictions of our approach are a peak-dip struc- 
ture of the spectral function, with a weakly dispersing 
peak at w ~ A and a peak-dip distance « A^. On the 
other hand we expect a broad incoherent peak which dis- 
perses like Ek/'^- In Fig-|3^ we present ARPES data for 
near optimally doped Bi2212 with Tc — 87K for mo- 
menta near a hot spotli£3. The intensity displays the 
predicted peak/dip/hump structure. A sharp peak is lo- 
cated at ~ 40meV, and the dip-as at SOmeV such that the 
peak-dip distance is 42meVlllZl. In the spin-fluctuation 
theory, the peak-dip distance is the energy of the INS 
resonance peak frequencyllZZHiZj. The neutron scattermg 
data on Bi2212 with nearly the same = 91K yieldE£^ 
As = 43meV, in excellent agreement with this prediction. 
Furthermore, with underdoping, the peak-dip energy dif- 
ference decreases and, within error bars, remains equal 
to Ag. This behavior is illustrated in Fig.Ba. 
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In Fig, 36 we present experimental results for the vari- 
ation of the peak and hump positions with the deviation 
from the Fermi surface. These show that the hump dis- 
perses with k — kp and eventually recovers the position 
of the broad maximum in the normal state. At the same 
time, the peak shows little dispersion, and does not move 
further in energy than A + As. Instead, the amplitude 
of the peak dies off as k moves away from kp. This be- 
havior is_aga.in fully consistent with the theoretical pre- 
dictionsEBB. 
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FIG. 35. The experimental peak-dip distance at various 
doping concentrations compared with As extracted from neu- 
tron measurements. Data from Ref.[30]. The theoretical re- 
sult is presented in Fig. 
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FIG. 36. The dependence of the experimental peak (flat 
curve) and hump (dispersing curve) positions on the devi- 
ation from the Fermi surface. The hump disperses with 
k — kp (dotted line) and eventually recovers the position of 
the broad maximum in the normal state, while the peak po- 
sitio n ch anges little with the deviation from kp. Data from 
Ref. 147. The theoretical result is presented in Fig. p^ . 



We regard the presence of the dip at A -f As, and the 
absence of the dispersion of the quasiparticle peak as two 
major "fingerprints" of strong spin-fluctuation scattering 
in the spectral density of cuprate superconductors. 



3. The density of states: 

The fermionic DOS N{llj) is proportional to the dy- 
namical conductance dl/dV through a superconductor- 
insulator- normal metal |(SIN) measured ai uj = eV where 
V is the applied voltagellHJ. The key prediction of our ap- 
proach is the occurrence of a dip in the DOS at an energy 
« As away from the peak at w = A. The drop in the 
DOS at 17t = A -f- As from Eq.|8^ can be understood in 
terms of SIN conductance as follows: when the applied 
voltage, V, equals flt/e an electron that tunnels from a 
normal metal can emit a spin excitation and fall to the 
bottom of the band, losing its group velocity. This loss 
leads to a sharp reduction of the current and produces a 
drop in dl/dV. This process is shown schematically in 
Fi^. 

SIN tunneling experimentSpiiave been performed on 



YBCO and Bi2212 materialsE£fl. Sinular results have 



been recently obtained by Davis et al.li^. At low and 
moderate frequencies, the SIN conductance displays a 
behavior which is generally expected in a d— wave super- 
conductor, i.e., it is linear in voltage for small voltages, 
and has a peak at eV p— A^-where A is the maximum 
value of the d— wave gapt^^t^ The value of A extracted 
from tunneling agrees well with the maximuw— walue of 
the gap extracted from ARPES measurementgi^'t^. At 



frequencies larger than A, the measured SIN conductance 
displays an extra dip-hump feature which become visible 
at around optiaial doping, and grows in amplitude with 



underdoping 
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At optimal doping, the distance be- 



tween the peak at A and the dip is around 40meV. This 
is consistent with As extracted from neutron measure- 
ments. 



4- SIS tunneling: 

The major prediction of the spin-fermion model for the 
SIS tunneling conductance, S{uj), is the emergence of a 
singularity at lo = 2 A -I- As. As mentioned above, this 
singularity is likely softened due to thermal excitations 
or non-magnetic scattering processes and transforms into 
a dip slightly below 2 A -I- As, and a hump at a frequency 
larger than 2 A -I- As. 
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data as an additional strong argument in favor of the 
magnetic scenario for superconductivity. 



5. Optical and Raman response: 




FIG. 37. SIS tunneling conductance normalized by A for 
Bi2212 materials ranging from overd oped (top curves) to un- 



derdoped (bottom curves) from Ref. 157, The peak dip dis- 
tance increases for increasing doping and saturates at around 
3A as expected in our theory. The corresponding theoretical 
result is presented in Fig. 

Recently, Zasadzinski et al. obtained both new data 
and carefully examined their previous SIS tunneling data 
for a set of|Bi2212 materials ranging from overdoped to 
underdopedt£d. Their data, presented in Fig.^ show 
that in addition to the peak at 2A, the SIS conductance 
displays the dip and the hump at larger frequencies. The 
distance between the peak and the dip (whicK-a: 
mately equals Ag in the spin fluctuation modelEIi 
close to 2A in overdoped Bi2212 materials, but goes down 
with underdoping. Near optimal doping, this distance is 
around 40meV. For an underdoped, Tc — 74K, mate- 
rial, the peak-dip distance is reduced to about 30meV. 
These results are in qualitative and quantitative agree- 
ment with ARPES and neutron scattering data, as well 
as with the theoretical estimates. The most important 
aspect is that with underdoping, the experimentally mea- 
sured peak-dip distance progressively shifts down from 
2A. This downturn deviation from 2A is a key feature 
of the spin-fluctuation mechanism. We regard the exper- 
imental verification of this feature in the SIS tunneling 
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I) is 



Theoretical considerations show that optical measure- 
ments are much better suited than Raman measure- 
ments tXL- search for the "fingerprints" of a magnetic 
scenariolllj. For the optical conductivity we predict a 
singular behavior at energies 2A -t- As, 4A, 2A -I- 2 As, 
which can be amplified if one considers the second deriva- 
tive of conductivity via W{uj) = -^{ijj'R,ea'^^{uj)). Ev- 
idence for strong coupling effects in the optical conduc- 
tivity in superconducting cuprates has been reported in 
Refs. p^ , |T60| , p^p^ 
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FIG. 38. A comparison between theoretical and e xper imen 



tal results for the optical condu ctivi ty (from Ref. 
experimental data are from Ref |l58[. 
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We first discuss the form of (Ti(w) 
pare the theoretical result for (Ji^r^^ 
perimental data by Puchkov et al. 
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Eig. ^ we com- 
113 with the ex- 



Ma for optimally doped 
YBa2Cu306+5 in the superconducting state. The pa- 
rameters are the same as in the normal state fits. As 
the theoretical formula does not include the contribu- 
tions from the nodes, the comparison is meaningful only 
for ix) > 2 A. We see that the frequency dependence of 
the conductivity at high frequencies agrees well with the 
data. The measured conductivity drops at about lOOmeV 
in rough agreement with 2A -I- As which for A « BOmeV 
and As ~ 40meV is also around lOOmeV. The good 
agreement between theory and experiment is also sup- 
portive of our argument that the momentum dependence 
of the fermionic dynamics becomes irrelevant at high fre- 
quencies, and fermions from all over the Fermi surface 
behave as if they were at hot spots. 

We next consider the singularities in the frequency de- 
pendence of the conductivity in more detail and com- 
pare the theoretical and experimental results for W{uj) = 

{ujKca^^ (uj)) . The theoretical result for W{uj) is pre- 
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sented in Fig.g6| The experimental result for W{lu) in 
YBCO is shown in Fig.^. We see that the theoretical 
and experimental plots of W{lu) look rather similar, and 
the relative intensities of the peaks are at least qualita- 
tively consistent with the theory. We identify (see ex- 
planations below) 2A + As with the deep minimum in 
W{oj). This identification, that is consistent with the 
analysis of ai{uj), yields 2A-t- ~ lOOmeV. Identifying 
the extra extrema in the experimental W{uj) with 4A 
and 2 A -I- 2As, respectively, we obtain 4A ~ ISOmeV, 
and 2A -I- 2As ISOmeV. We see that three sets of data 
are consistent with each other and yield A ~ 30meV and 
As ~ 40 — 45meV. The value of in good agree- 
ment with tunneling measurements li^, and A^ agrees 
well with the resonaHce frequency extracted from neu- 
tron measurementst£2l. Indeed, the analysis of a second 
derivative of a measured quantity is a very subtle proce- 
dure. The good agreement between the theory and ex- 
periment is promising but has to be verified in furtbej; on-, 
perimental studies. Still, theoretical calculations |123'Ef2l 
clearly demonstrate the presence and observability of 
these "higher harmonics" of the optical response at 4 A 
and 2A + 2A,. 
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FIG. 39. 
mental results for W{lu) — -^j^lujllea^ ' (uj 
The theoretical result is presented in Fig. M. The position 
of the deep minimum agrees well with 2 A -|- As. The extrema 
at higher frequencies are consistent with 4A and 2 (A -|- As) 
predicted by the theory. 



Finally, we comment on the position of tl 



A, 

peak and compare the res«]ts of Abanov et alliii^lilii with 
those by Carbotte et di^. Theoretically, at T = and 
in clean limit, the maximum and minimum in W{uj) are 
located at the same frequency. At a finite T, howevpet, 
they quickly move apart (see Fig ^) . Carbotte et ali£^ 
focused on the maximum in W(uj) and argued that it is 
located at A-l-A^ instead of 2 A -I- A^. We see from Fig^ 
that the maximum in W{uj) shifts to a lower frequency 
with increasing temperature and over some T range is lo- 
cated close to A -I- As. On the other hand, the minimum 
in W{uj) moves very little with increasing T and virtually 



remains at the same frequency as at T = 0. This result 
suggests that the minimum in W{uj) is a more reliable 
feature for comparisons with experiments. This conclu- 
sion is in agreement witii, recent conductivity data on 
optimally doped Bi2212ll23. W{uj) extracted from these 
data shows a strong downturn variation of the maximum 
in W{uj) with increasing temperature, but the minimum 
in W{uj) is located at around llOmeV for all tempera- 
tures. 



D. Experimental facts that we cannot yet describe 

There are several experimental results that we do not 
understand. EiES|l_are the results by Ando, Boebinger and 
collaboratorsEiirEij on the behavior of the Lanthanum 
and Bismuthate based superconductors in magnetic fields 
sufficiently strong to (almost) destroy superconductivity. 
For doping levels close to the optimal one, they found 
that the resistivity at low temperatures continues to be 
linear in T with the same slope seen at higher tempera- 
tures. If the assumption that the magnetic field destroys 
superconductivity but otherwise does not affect the sys- 
tem properties is correct, this result poses a problem for 
the spin-fluctuation model as the latter yields a linear in 
T resistivity over a wide range of temperatures, but only 
for T larger than a fraction of Wsf. To account for these 
data one might have to invoke some kind of quantum- 
critical physics associated with the opening of the pseu- 
dogap (see below). 

Another experiment that is not yet understood is the 
measurement of the Hall angle, 9h = Pxy/ Pxxr-^ii^ch 
shows an incredibly simple behavior, cot Oh oc T^po"^! pnH 



also displays a particular frequency bchaviorEHj'ElH 
The orbital magnetoresistance Ap/p also behaves in quite 
an unusual way, violating Kohler's rule, according to 
which Ap/p is a function of H^/p^, independent of T, 
where H is the applied magnetic field. Some of this 



lenomenological 
however prob- 



physics is already captured in the semi- 
calculations by Stojkovich and Pinctcia; 
lems remain. In the description based on the spin- 
fermion model the technical problem not yet solved is 
how to include in a controlled way vertex corrections 
which are not small; in one of the vertices for the Hall 
conductivity the momentum transfer is small. Some 
progress with these calculations, have been recently made 
by Katami and coUaboratorsEHI. Another explanation of 
the Hall data, has recently been proposed by Abrahams 
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and Varmac 

Yet another unanswered question, already noted 
above, is the origin of a large (almost lOOmeV), frequency 
and temperature independent contribution to the self- 
energy that one has to invoke in order to fit conductivity 
and ARPES data. It could, in principle, be due to in- 
elastic scattering by impurities, but its very large value 
makes this explanation problematic. 

Electronic Raman scattering reveals further puzzling 
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behavior: in all geometries one observes a frequency in- 
dependent behavior over a very large energy scale, fre- 
quently referred to as the positive background. More- 
over, the overall size o f the background is very different 
in different geometriealI3llIj. 

There are also uncertainties associated with reconcil- 
ing the incommensurability of tbe-pagnetic response in 
the normal state of 214 matcrialalla with the commensu- 
rate peaks required to obtain a consistent explanation of 
-'^^O and ^^Cu NMR experiments, but these are not likely 
to pose fur^flai jnental problems to the spin-fluctuation 
approach 
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Finally, the claim of universality of the low-energy 
behavior relies heavily on the existence of a quantum 
critical point at which the antiferromagnetic correlation 
length diverges. In real materials there are indications 
that the transition to antiferromagnetism is actually of 
first order. In this situation, the theory we described 
is valid only if there still exists a substantial region in 
parameter space where the system is critical before it 
changes its behavior discontinuously. NMR and neutron 
scattering experiments on optimally doped cuprates seem 
to support such behavior. Another reason for concern is 
the role of disorder and inhomogeneities. Despite enor- 
mous progress in sample fabrication, cuprates often tend 
to be very heterogeneous materials. It has been estab- 
lished in several cases that these aspects are actually 
intrinsic, forcing one to include effects due to inhotan 
geneities and disorder into the theoretical descriptionc^^ 



E. Phase diagram 

In this section we discuss in detail the experimental 
phase diagram of cuprate superconductors and comment 
on the origin of the pseudogap behavior found for small 
charge carrier concentrations. 

From a general perspective, the key to understanding 
of cuprate superconductors is identifying the nature of 
the protected behavior of the novel states of matter en- 
countered in the insulating, conducting, and supercon- 
ducting states as one varies doping and temperature, 
including the possible existence of one or more quan- 
tum critical points. Consider first the YBa2Cu307_5 
system oii— which the generic phase diagram of Fig.|l| 
was basedll^. A somewhat similar diagram based on 
transport mcaaurements was independently proposed by 
Hwang et alM3, while one based on specific heat and sus- 
ceptibility rneasurements has been proposed recently by 
Tallon et al.t^a. As discussed in the Introduction, in ad- 
dition to the Tc line, there are two crossover or phase 
transition lines in Figjl]. The upper line T — Tcr is 
defined experimentally by a maximum in the tempera- 
ture dependent uniform magnetic susceptibility, xo- It 
has been further characterizedll23 as the temperature at 
which the antiferromagnetic correlation length ^ is of the 
order Cu-Cu lattice spacing (Barzykin and Pines used a 



criterion ^(Tcr) = 2a). The lower line T = T* may be 
defined experimentally as the temperature at which the 
product of the copper spin-lattice relaxation time,^^Ti 
and the temperature, T, reaches its minimum value. In 
the Bi2Sr2CaCu408 counterparts of the YBa2Cu307_5 
system, it corresponds to the temperature at which the 
leading edge gap found in ARPES experiments for quasi- 
particles near (tt, 0) becomes fully open, effectively gap- 
ping that portion of the quasiparticle Fermi surface. To 
a first approximation, on making use of the experimen- 
tal results for optimally and underdoped YBa2Cu307_5 
materials one finds that 
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The superconducting in Fig. |l] is obtained using the 
empirical relation t2B 



0.16±0.11, 



Tc 
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where x is the doping level, and T^^^ is the maximal 
transition temperature for a given class of materials. The 
location of Tcr can well be fitted by another empirical re- 
lation. 
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where Xcr ~ 0.19. Similar expressions are found for the 
La2_2;Sr3;Cu04 and Bi2Sr2CaCu408 materials. This ex- 
pression for Tcr is, in both its magnitude and doping de- 
pendence, close to the pseudogap-temperature obtained 
by Loram and his coUaboratoraiHl from an analysis of 
specific heat experiments. A remarkable result of this 
purely phenomenological analysis is that the crossover 
temperature Tcr extrapolates at zero doping to the known 
value of the antiferromagnetic super-exchange interaction 
J. 

The fit to Tcr by Eq. 104 raises the issue of whether 



T* and Tcr are independent of Tc and would extrapolate 
to the origin at a doping level x = Xcr if superconduc- 
tivity was absent. The system then would have an ad- 
ditional quantum critical point at a; = Xcr with a new 
kind of ordered state for x < x^- This issue is cur- 
rently open and is a subject of active research. Support 
for a phase diagram with an additional quantum critical 
point at X = Xcr comes from the work of Loram , Tallon, 
and their collaboratorsllHriHa, who have proposed such 
behavior based on a detailed analysis of their specific 
heat experiments on underdoped and overdoped systems. 

Iknner, Panagopoulos and others 



Moreover, as Loram. 

have emphasized E23 t£j, in the superconducting state 
of the low-doping side of T^ one has "weak" supercon- 
ducting behavior, with a superfluid density ps decreasing 
with decreasing doping, while on the high doping side 
one has a "conventional" superconductivity, and a value 
of Ps that is nearly independent of the doping concen- 
tration. Further support for the idea of an additional 
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quantum-critical point comes from the well established 
fact that optimally doped cuprates are the ones for which 
the extension of the linear resistivity to T = yields very 
small residua l resistivity, and from the experiments of 
Refs. 214,215 which, we recall, show that in the absence 



of superconductivity the linear temperature dependence 
of the resistivity extends to lower T indicating that at 
some doping the resistivity can be linear down to T = 0. 
As Laughlin et a/.c22l have emphasized, the presence of 
a quantum critical point with a large domain of influ- 
ence, together with superconductivity, serves to conceal 
the nature of the non-superconducting ground states on 
either side of the quantum critical point. One might hope 
that ARPES experiments near optimal doping would dis- 
tinguish between a quantum critical behavior with quan- 
tum critical point at around optimal doping and a spin 
fermion scenario with antiferromagnetic quantum criti- 
cal point at considerably smaller doping concejjtration. 



However, a recent analysis of Haslinger et aZ.liiil showed 
that fits to current experiments with either model is pos- 
sible and requires in both cases the introduction of a large 
temperature and frequency independent scattering rate, 
as noted earlier. 

The variety of experimental results for the pseudogap 
allows one to understand it phenomenologically, without 
invoking a particular microscopic mechanism. First, as 
T* and Tcr scale with each other, it is natural to at- 
tribute both T* and Tcr to different aspects of the same 
physical phenomenon which begins at Tcr and gains full 
strength at T* . This idea is fully consistent with NMR 
data which show the onset of changes in TiT at T^r, which 
eventually give rise to a sign change of the temperature 
derivative of ^'^TiT at T* . Second, ARPES data on the 
leading edge gap clearly demonstrate that the pseudogap 
physics below T* is associated with the redistribution of 
the spectral weight for hot quasiparticles; quasiparticles 
near the nodes are almost unaffected by the development 
of the pseudogap. In the ARPES literature, this effect is 
described as a progressive development of the arcs of the 
Fermi surface centered around nodal points. The evo- 
lution of the full Fermi surface into the arcs begins at 
around T*, and at Tc the whole Fermi surface becomes 
gapped. The "gapping" of hot fermions obviously affects 
NMR relaxation rates dominated by momenta near Q 
(such as ^^Cu TiT) as a spectral weight transfer would 
lead to a reduction of a deea#-|rate of a spin fluctua- 
tion into a paxticle-hole paiiE3^E§ — NMR experiments by 
Curro et a/.E12l and Haase et alM3 show that this is in- 
deed the case. The gapping of hot quasiparticles should 
also lead to a temperature-depeadent reduction in the 
uniform magnetic susceptibilitylll3. 

The phenomenological description is of course not 
enough as it leaves open the key question, namely what 
causes the spectral weight transfer for hot quasiparticles. 
We now discuss how the experimental phase diagram fits 
into the spin-fluctuation scenario. 

First of all a general phase diagram based on spin fluc- 
tuation approach should distinguish between weakly and 



strongly antiferromagnetic materials. Weakly antiferro- 
magnetic materials are those to the right of Tcr for which 
the dimensionless coupling constant A is smaller than 
unity, which corresponds to a correlation length smaller 
than a few lattice constants. For these materials, the 
normal state is a renormalized Fermi liquid, the nearly 
antiferromagnetic Fermi liquid, and Tc signals a transi- 
tion to a BCS-like superconducting state with a (1^2 _y2 
order parameter. For strongly antiferromagnetic materi- 
als, on the other hand, the dimensionless coupling con- 
stant A > 1 and ^ > 2a. In this situation, the normal 
state behavior deviates from a Fermi liquid already at 
comparatively small w and T although at the lowest fre- 
quencies the system still would display a Fermi liquid be- 
havior if indeed one could extend the normal state down 
to r = 0. For these systems, we also know that the pair- 
ing instability temperature Tcr increases with decreasing 
doping and for large enough coupling saturates at a value 
comparable to the magnetic J. Applying this to the ex- 
perimental phase diagram, we see that optimally doped 
materials are at the borderline between being weakly and 
strongly antiferromagnetic: on the one hand the dimen- 
sionless coupling is already not small, on the other hand, 
the pseudogap phase extends at best over a T range which 
is only a fraction of Tcr • 

As wc have noted above, these results of the spin- 
fermion model make it a strong candidate for the micro- 
scopic description of the pseudogap phase: Tcr saturates 
Tcr saturates at a finite value at the magnetic transi- 
tion; for A > 1, the pairing involves non- Fermi liquid 
fermions; at T = there are two distinct energy scales 
in the problem, a fermionic gap A oc Tcr and a bosonic 
gap (X TcrA~^ <C Tcr- The central issue is whether the 
pairing of incoherent fermions only gradually changes the 
fermionic self-energy, or whether it creates a feedback 
on fermions which immediately gives rise to a coherent 
quasiparticle behavior at the lowest frequencies, as hap- 
pens in dirty superconductors where S(ti-') = i"f in the 
normal state transforms below Tc into a mass renormal- 
ization at the smallest lu 



I](a;) = i"f- 



A2)i/2 
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If the feedback is gradual, then the pairing creates bound 
states of incoherent fermions with S = 0. In this situa- 
tion, there is a reduction in the density of states below 
Tcr, but a full superconducting gap does not develop un- 
til a smaller temperature, Tc- A sirupJe toy model which 
describes this physics would be onec3 in which fermions 
in the normal state display a quantum critical behav- 
ior with S(a') « S(ti;) = (icjw)^/^, and pairing creates a 
nonzero flat pairing vertex $ but does not affect S(w). 
In this situation, the fermionic propagator acquires a gap 
at a finite but imaginary frequency: 



Gk(w) cx 



iuj — 



(106) 



where 
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The spectral function at k = \ip and the fermionic DOS 
both have broad maxima at w = A = ^^/cj, but the 
spectral weight is finite at any finite w, although reduced 
at low frequencies. Within this model, the transition to 
the true superconducting st ate c an be understood as a 
rotation of the pole in Eq. (|lO(j ) from the imaginary to 
the real frequency axis such that at the lowest frequen- 
cies the iuj term becomes purely real. The frequency up 
to which N{u}) = then give an estimate for the actual 

It is unclear to what extent the results of this toy model 
reflect the physics of the spin-fermion model below Tcr- 
Without elaborating on this subject of current research 
we mention that as long as an Eliashberg approach is 
justified, phase fluctuations of the superconducting or- 
der parameter cannot substantially reduce Tc compared 
to Tcr- Behavior different from that in dirty supercon- 
ductors could emerge only if longitudinal fluctuations of 
the superconducting order parameter are soft and able 
to destroy the superconducting coherence s^^c ^ Tcr. 
We have both numericalE23 and analytical c33 evidence 
that such degeneracy does exist in the limit A = oo. Still, 
this subject is far from being fully understood and clearly 
requires further study. 

Another subtle issue is whether the spin-fermion model 
displays quantum critical behavior at a; = Xcr- Phys- 
ically, this would imply that the pairing of incoherent 
fermions at Td- and the pairing of coherent fermions at Tc 
are uncorrelated phenomena - the first gives rise only to 
the pseudogap, while the latter yields BCS-type super- 
conductivity. Since the pairing of incoherent fermions 
is not a perturbative phenomenon aiid-|requires the in- 
teraction to exceed a threshold valueE£3, the pairing of 
incoherent fermions involves only quasiparticles in some 
finite region around a hot spot and would form a dome 
on top of a magnetic quantum critical point and vanish 
at a finite x. It is not clear how well one can separate 
coherent and incoherent pairings. 

Another issue related to the possible explanation of 
the pseudogap within the spin-fermion model is whether 
one can smoothly interpolate between the limit g <^ W 
(W is the fermionic bandwidth), where one can perform 
calculations in a controlled fashion, and g ^ W where 
Mott physics become relevant. In essence the issue is 
whether or not there is a qualitative difference between 
limits in which the effective interaction is either much 
larger or much smaller than the fermionic bandwidth. 
In the latter case, it is appealing to conjecture that the 
pseudogap is associated with the fact that it is difficult 
for hot quasiparticles to be both itinerant and localized. 
Tcr then would mark the onset of insulating behavior as- 
sociated with such localization, and the pseudogap phase 
would represent a kind of partial Mott insulator. A sign 
that these two limits may describe some aspects of the 
basic physics similarly is that, as shown earlier, Tcr scales 
with g ioT g < W, but crosses over to /g oc J for 



g > W. At the same time, the limit ^ > is probably 
more rich than the small g limit as the physics associ- 
ated with the localization in the Mott insulator is not 
included into our analysis. We speculate that due to this 
Mott physics, bound singlet pairs of fermions that emerge 
below Tcr could order for example in coiyapnar fashion as 
suggested by Sachdev and collaboratorg222l who arrived at 
a columnar phase by studying weakly doped Mott antifer- 
romagnets. This ordering in turn would imply that Tcr is 
a true phase transition line below which Z4 symmetry is 
broken. Columnar ordering also opens a link between our 
approach and the approaches which depart from Mott in- 
sulator at half-filling. In particular, columnar ordering-of 
bound electron pairs naturally leads to stripe physicsE2£ . 
An alternative p,©savbility is that singlet pairs remain spa- 



tially disordere(£l3^cH3. In any event, the role of localiza- 



tion effects certainly increases as the system approaches 
half-filling. Whether they remain strong near optimal 
doping in the normal state is a subject of debate, but 
still, strong localization effects should reduce the number 
of low energy carriers and therefore change the volume 
of the Fermi surface or increase their mass. ARPES ex- 
periments on the other hand indicate that in the normal 
phase, the Fermi surface is large and obeys Luttinger's 
theorem without dramatic mass renormalizations. We 
therefore believe that near optimal doping localization 
effects are at best moderate. 

We conclude this discussion of the phase diagram by 
mentioning two alternative scenarios for the pseudogap 
and anomalous normal state propepties. The first sce- 
nario, pioMered by P.W^ AwjiersonEia, X.-G. Wen and 
P. A. Leeuij and othersE2°rE2Zl and. later modified chiefiy 
by M.P.A. Fisher and T. SenthilE23, assumes spin-charge 
separation at half-filling and explains the whole phase 
diagram as a result of weak doping of a Mott insula- 
tor. A second scenario, on the contrary, assumes that 
one can understand the phase-diagram within mean-field 
theory; as Chakravarty et al£3 have proposed the pseu- 
dogap might then be a new protected state of matter that 
is the result of the breaking of a hidden symmetry. 



VII. CONCLUSIONS 



In this Chapter we have demonstrated that supercon- 
ducting pairing mediated by the exchange of spin fiuc- 
tuations is a viable alternative to conventional phonon- 
mediated pairing. We discussed in detail the normal 
state properties, the pairing instability and the super- 
conducting behavior of a material near an antiferromag- 
netic instability, when the dominant interaction between 
quasiparticles is of electronic origin and, at energies much 
smaller than the fermionic bandwidth, can be viewed as 
being due to the emission and absorption of a collective, 
soft spin degree of freedom. We argued that the spin- 
fluctuation exchange yields an attraction in the (1^2 _y2 
channel in agreement with what nearly all researchers 
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now believe is the pairing symmetry in the cuprates. 

We demonstrated that the physics is quahtatively dif- 
ferent depending on whether or not Fermi surface geom- 
etry allows a process in which a collective mode decays 
into a particle and a hole. For this to be possible, the 
Fermi surface should contain hot spots. Wc focused on 
the case in which the Fermi surface does contain hot spots 
(as the photoemission experiments in cuprates indicate) 
and showed that spin fluctuations are then overdamped 
and that their diffusive dynamics should be analyzed in 
a consistent manner with the low energy dynamics of the 
fermions. We further argued that contrary to naive ex- 
pectations, this case is better for d— wave pairing than 
one in which spin fluctuations are propagating, magnon- 
like quasiparticles. 

We showed that the low-energy theory for fermions in- 
teracting with overdamped collective spin excitations is 
universal, independent of the details of the underlying 
lattice Hamiltonian and is characterized by only two in- 
put parameters: the dimensionless coupling constant. A, 
and an overall energy scale, w, proportional to the effec- 
tive spin-fermion interaction In so doing, we have de- 
veloped the microscopic justification for the NAFL. The 
coupling constant. A, scales with the magnetic correlation 
length, so that close enough to a magnetic transition, the 
system falls into a strong coupling regime. A universal 
description is valid if g is smaller than the fermionic band- 
width W (that would correspond to a weak coupling limit 
if the system was far away from a magnetic instability). 
In the opposite limit which we did not discuss in detail, 
lattice effects become important, and the universality is 
lost. 

At sufficiently low temperatures and energies, the 
nearly antiferromagnetic Fermi liquid (NAFL) is a Fermi 
liquid quantum protectorate, according to Landau's 
criterion-that one can obtain a one-to-one correspon- 
dence between the low-lying states of a Fermi gas and 
the Fermi liquid, as though the particle interaction was 
turned on adiabaticallyll^. However, because of the 
closeness to an antiferromagnetic instability, it is an un- 
conventional Fermi liquid, in that the characteristic en- 
ergy above which this description is no longer valid is 
not the Fermi energy, but is the much much lower spin- 
fiuctuation energy Wgf = (D/(4A^) that is typically two 
orders of magnitude smaller than the Fermi energy. For 
energies (or temperatures) between Ugf and the Fermi en- 
ergy, we have seen using the spin-fermion model that the 
system behavior is again universal and depends only on 
a very limited number of phenomenological parameters. 
So in this sense, the behavior of an NAFL at energies and 
temperatures above cogf is also protected, and one finds, 
in the NAFL, two distinct protected states of matter, 
depending on the energy or temperature one encounters. 

We compared in detail the spin fluctuation approach 
with the Eliashberg approach to phonon superconductors 
and showed that despite the absence of the small elec- 
tron to ionic mass ratio that justified Eliashberg theory 
for phonons, an Eliashberg-type approach to the spin- 



fermion model is still justified, but only at strong cou- 
pling A > 1. 

We showed that at large A, there are two distinct 
energy scales for the normal state problem: Q — Ag, 
A — 0(1), which at g ^ is of the order of (even 
though numerically smaller by about an order of mag- 
nitude than) the Fermi energy, and a much smaller 
LOsi = i'/(4A^) <C UJ- Conventional Fermi liquid behavior 
with S" (X uj'^ + (ttT)^ and an almost temperature inde- 
pendent static spin susceptibility exists only at frequen- 
cies and temperatures smaller that Wgf- At frequencies 
between tUgi and ZU, the system crosses over into a regime 
in which spin susceptibility is diffusive and E" oc to". We 
found that a — ^, while the behavior in the crossover 
region from Fermi liquid to non Fermi liquid behavior 
resembles a linear in frequency dependence. In the dif- 
fusive regime, both spin and fermionic propagators be- 
come independent of A, i.e., the system displays magnetic 
quantum-critical behavior. Finally, at w > U, the self- 
energy gets smaller than the bare w, although S" oc cj^/^ 
still holds. 

We next argued that the existence of the two dis- 
tinct energy scales has a strong impact on the pairing 
problem: the pairing instability temperature Tcr is pre- 
dominantly determined by incoherent, non- Fermi liquid 
fermions with energies between Wgf and uj. We then con- 
sidered the superconducting state, found the pairing gap 
A oc Tcr (2A/Tcr 4 for A > 1), and argued that one 
of the signatures of the spin-fiuctuation mechanism is 
the change of the low-frequency spin dynamics from re- 
laxational to propagating due to the feedback from the 
pairing. This give rise to the emergence of the resonance 
peak in the spin propagator at a frequency As which at 
large A scales as As l.STcr/A. For smaller A, the reso- 
nance frequency goes up but always remain smaller than 
2 A. This restriction (not found for phonon superconduc- 
tors) is a consequence of feedback between fermionic and 
bosonic dynamics. 

We discussed how the emergence of the resonance 
mode affects fermionic properties and identified sev- 
eral "fingerprints" of the strong coupling spin-fluctuation 
scenario. These include, but arc not limited to, the 
peak/dip/hump in the spectral function, the near- 
absence of the dispersion of the quasiparticle peak, the 
peak/dip features in the SIN and SIS tunneling conduc- 
tances, and in the second derivative of the optical con- 
ductivity. 

Finally, we compared the theoretical results with the 
normal state data for near optimally-doped cuprates, and 
found that a large number of experimental results can be 
understood within the spin-fluctuation scenario. In par- 
ticular, we argued that the tunneling, photoemission and 
optical data for the superconducting state display the ex- 
pected "fingerprints" of the spin-fluctuation mechanism 
for superconductivity in the cuprates. 

The spin-fluctuation theory should be equally appli- 
cable to quasi two dimensional organic superconductors 
or compounds of the family CeXIns with X^Co, Rh, 
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Ir, if the Fermi surface of these systems possesses hot 
spots. A detailed experimental investigation of spectro- 
scopic properties of these systems in the superconducting 
state will provide important clues on whether they are 
indeed magnetically mediated superconductors. 

Despite the successful identification of the "finger- 
prints" of spin mediated pairing in the cuprates, we listed 
in Section 6.4 some properties of these materials which 
we have not explained within the spin-fluctuation sce- 
nario. Perhaps the most important remaining question 
is to what extent the spin-fluctuation theory explains the 
pseudogap physics in underdoped cuprates. The exis- 
tence of the two distinct energy scales in the spin-fermion 
model and the presence of fermionic incoherence in the 
normal state make this model a good candidate model 
for the pseudogap. The plausible argument here is that 
while the pairing of the incoherent fermions will create 
spin singlets, fermions will still behave incoherently and 
hence not superconduct until a smaller is reached 
when feedback from the pairing restores fermionic co- 
herence. Particularly relevant here is whether or not 
there exists an additional critical point on the phase di- 
agram at around x = 0.19. This critical point emerges 
in the spin-fluctuation theory if the pairing of incoherent 
fermions from frequencies above ujsf and the pairing of 
coherent fermions from frequencies below ujsf are sepa- 
rate phenomena (the first then definitely gives rise to a 
pseudogap physics while the second yields a true super- 
conductivity). The pairing of incoherent fermions is a 
threshold phenomenon and thus occurs at T^-cix) ■ We 
still however need to understand whether incoherent pair- 
ing and coherent pairing can be totally separated from 
each other. 

Another unresolved issue is to what extent the fact 
that the parent compounds of the cuprates (La2Cu04 or 
YBa2Cu306) are Mott insulators and nearest-neighbor 
Heisenberg antiferromagnets affects the behavior of 
doped materials. For the spin-fermion model, this ques- 
tion could be reformulated as whether the limit g <C W, 
in which we can separate low and high energies and per- 
form controlled calculations, and the limit g ^ W, are 
only quantitatively different or are qualitatively differ- 
ent. An encouraging sign that the two limits may not be 
very different is the result that Tcr smoothly interpolates 
between 0{g) for g < to olw'^/g) for 5 > PF, in 
the latter case Tcr becomes of order of a magnetic ex- 
change integral J. On the other hand, the localization 
effects which accompany a transition to a Mott insulator 
are not included in the theory presented here. In partic- 
ular, we cannot predict what would happen with singlet 
pairs below Tcr, if it indeed is the onset of the pseudo- 
gap - whether they remain disordered-or form columnar 
stripes as Sachdev and collaboratorsc23 suggested. 

The enduring presence and richness of unsolved prob- 
lems in the field of unconventional superconductivity 
makes us optimistic that we will continue to see unex- 
pected experimental observations and new, creative theo- 
retical concepts in this field of research. We have enjoyed 



over a number of years our respective collaborations on 
the spin fluctuation approach because of its clarity and 
ability to make falsifiable predictions for experiment. We 
hope that the readers of this Chapter will share our ex- 
citement for this approach to magnetically mediated su- 
perconductivity. 
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